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Abstract 

We consider population models in which the individuals reproduce, die and also 
migrate in space. The population size scales according to some parameter N ^ which can 
have different interpretations depending on the context. Each individual is assigned 
a mass of and the total mass in the system is called population density. The 
dynamics has an intrinsic density regulation mechanism that drives the population 
density towards an equilibrium. We show that under a timescale separation between the 
slow migration mechanism and the fast density regulation mechanism, the population 
dynamics converges to a Fleming- Viot process as the scaling parameter N approaches 
oo. We first prove this result for a basic model in which the birth and death rates can 
only depend on the population density. In this case we obtain a neutral Fleming- Viot 
process. We then extend this model by including position-dependence in the birth and 
death rates, as well as, offspring dispersal and immigration mechanisms. We show how 
these extensions add mutation and selection to the limiting Fleming- Viot process. All 
the results are proved in a multi-type setting, where there are q types of individuals 
interacting with each other. We illustrate the usefulness of our convergence result by 
discussing applications in population genetics and cell biology. 
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1 Introduction 

Density-dependent models are well-known in population biology. In these models, the birth 
and death rates of individuals may depend on the density of the population, where the term 
density refers to the population size under a suitably chosen normalization. Many models in 
ecology, epidemiology and immunology can be suitably described by such models (see Thieme 
[35]). Considering the molecules of chemical species as the individuals in a population, we 
can also view a chemical reaction network as a density-dependent population model. 

Density-dependent models are appealing because one can easily account for interactions 
among individuals by appropriately specifying the birth and death rates as functions of 
the population density. For competitive interactions, as in the Lotka-Volterra model (see 
[28, 38, 30]), the death rate increases with the population density while for cooperative 
interactions, as in the Allee model (see [1]), the birth rate increases with the population 
density. Density dependent models are good candidates for modeling natural populations 
that cannot grow indefinitely due to the limited availability of certain vital resources or due 
to severe competition at large population sizes. By having the death rate dominate the birth 
rate at large densities, one can ensure that the population density does not go beyond a 
certain threshold. 

For a population having q types of individuals, the population density is a g dimensional 
vector whose i-th component gives the density of the population of the i-th type. For such a 
multi-type population, a density-dependent model can be written in the deterministic setting 
as a system of q ordinary differential equations. If all the trajectories of this system stay 
within a compact set at all times, then we say that the population dynamics has a density 
regulation mechanism. Such a mechanism is called equilibrating if all the trajectories reach 
a fixed point for this system as time goes to infinity. In such a situation, this fixed point is 
called the equilibrium population density. 

In this paper, we will consider population models in which the individuals live in a geo- 
graphical region E, that is a compact metric space. Even though we have a spatial structure, 
for us the population density will always denote the population size divided by a normal- 
ization parameter N. In other words, our notion of population density is global in the sense 
that it carries no information about the distribution of individuals in E. This is unlike other 
models of spatial populations where the population density is a spatially varying function 
specifying the local concentration of individuals at each location. The normalization param- 
eter N will be a large positive integer which can have various interpretations depending on 
the context. In ecological models, N can be taken to be the carrying capacity of a habitat, 
which is the maximum number of individuals that the habitat can support with its resources. 
In epidemic models, is usually the total population size, while in chemical reaction net- 
works, measures the volume of the system. In each of these cases, the population size at 
any time is of order A^. 
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For the moment assume that all the individuals have the same type and that they repro- 
duce, die and also migrate in E. At the time of birth, the offspring gets the same location as 
its parent. The population consists of approximately individuals of a mass of 1/iV each. 
The population density at any time is just the total mass of the individuals that are alive. 
Suppose that the birth and death rates of the individuals depend on the population density 
in such a way that they induce a density regulation mechanism which is equilibrating. We 
also assume that the migration mechanism operates at a timescale that is times slower 
than the density regulation mechanism. In such a setting, we can view the dynamics of the 
empirical measure of the population as a measure-valued Markov process parameterized by 
A^. Our goal is to understand how this family of Markov processes behaves as A^ — > oo. The 
population dynamics has two timescales separated by A^. If we observe the process at the 
fast timescale, then the effect of migration vanishes in the limit and it is uninteresting to 
consider the population with a spatial structure in this case. Therefore we will observe the 
dynamics at the slow timescale and examine its behaviour in the infinite population limit. 
Since the density regulation mechanism is fast, it will have enough time to re-equilibrate 
the population density between any two events at the slow timescale. Hence in the limit 
A^ — oo, we would expect the population density to remain equilibrated at all times. We 
will show that it is indeed the case. However our main task is to understand the dynamics of 
the spatial distribution of the population in the infinite population limit. We will prove that 
in the limit, the spatial distribution of the population evolves according to a Fleming- Viot 
process which takes values in the space of probability measures over E. This process was 
introduced in the context of population genetics by Fleming and Viot [17] in 1979 and it has 
been very well-studied since then. An excellent survey of Fleming- Viot processes is given by 
Ethier and Kurtz [13]. The model we just described will be called the basic model. In this 
model, the birth and death rates of an individual were density dependent, but independent of 
the location of the individual. It shall be seen later that the limiting Fleming- Viot process 
in this case is neutral, in the parlance of population genetics. If we add small position- 
dependent terms to the birth and death rates, then in the limit we obtain a Fleming- Viot 
process with genetic selection. Perhaps unsurprisingly, altering the birth rate this way leads 
to fecundity selection, while altering the death rate leads to viability selection, in the lim- 
iting Fleming- Viot process. We also consider extensions of the basic model by allowing for 
offspring dispersal (offspring is born away from the parent) or immigration. Such extensions 
add extra mutations to the limiting Fleming- Viot process. 

The results mentioned in the previous paragraph are proved in a multi-type setting. 
The population has q types of individuals and each type of individual can give birth to an 
individual of each type. All the individuals are migrating in E according to a type-dependent 
mechanism. Now we can view the joint dynamics of the empirical measures of the q sub- 
populations as a Markov process parameterized by A^. We make similar assumptions on 
the dynamics as before. Again in the limit A^ — )■ oo, the population density (which is now 
a g-dimensional vector) stays at an equilibrium at all times. Assuming the irreducibility 
of an underlying interaction matrix, we show that in the limit all the q sub-populations 
become spatially inseparable. This means that on any patch of E, either there is no mass 
present or there is mass of each type present in a proportion determined by the equilibrium 
density. Moreover the spatial distribution of each of the q sub-populations evolves according 
to a single Fleming- Viot process. This Fleming- Viot process can be seen as describing the 
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limiting dynamics of a mixed population, formed by taking a suitable density-dependent 
convex combination of the q sub-populations. 

In ecological models, the individuals need resources to survive and reproduce. Normally 
in spatial population models, resources are assumed to have a fixed distribution in space. 
As individuals move, they find the unexploited resources and compete for them locally with 
other individuals present in their neighbourhood. Such a model is different from the models 
we consider in two ways. Firstly, due to the local nature of the interactions, the density 
is locally regulated rather than globally regulated as in our models. Secondly, since the 
discovery of resources is tied to the movement of individuals, it is reasonable to assume that 
both migration and birth-death mechanisms operate at the same timescale. For such spatial 
models, Oelschlager [31] has shown in a multi-type setting that the dynamics converges in the 
infinite population limit to a system of reaction-diffusion partial differential equations. Such 
equations are in widespread use in biology (see Fife [16]). We now discuss the conditions 
under which our models can be useful. Consider a situation where the resource is not fixed 
but rapidly mixing in the whole space. This resource is shared by all the individuals in the 
population. An individual may deplete the resource locally but its effect is felt globally due 
to the rapid mixing. This gives rise to global density dependence in a spatial population. 
If the individuals move very slowly in comparison to their resource consumption mechanism 
(which is linked to their birth and death mechanisms), then we have a situation in which 
our models can be used. 

This paper is motivated by our earlier work [19] in which we study the phenomenon of 
cell polarity using a model considered here. Cell polarity refers to the clustering of molecules 
on the cell membrane. This clustering is essential to trigger various other cellular processes, 
such as bud formation [1] or immune response [39]. Therefore understanding how cells 
establish and maintain polarity is of vital importance. In [2], Altschuler et. al. devised a 
mathematical model for this phenomenon, by abstracting the mechanisms that are commonly 
found in cells exhibiting polarity. Their model has a fixed number of molecules that can either 
reside on the membrane or in the cytosol. These molecules move slowly on the membrane 
but diffuse rapidly in the cytosol. The dynamics has a positive feedback mechanism which 
allows a membrane molecule to pull a cytosol molecule to its location on the membrane. This 
mechanism is like a birth process in which a membrane molecule gives birth by exploiting 
the common resource (cytosol molecules) shared by all the membrane molecules. Since the 
migration of membrane molecules is slow and the mixing of the resource is fast, this model 
can be viewed as a model described in this paper (see Section 3.2 for details). Therefore 
the results in this paper are applicable and we obtain a Fleming- Viot process in the infinite 
population limit. In [19] we prove this convergence^ and use the limiting process to answer 
some interesting questions about the onset and structure of cell polarity. The model studied 
in [19] is rather simplistic as all the molecules are assumed to be identical. Most cells 
that exhibit polarity have molecules of many different types participating in the feedback 
mechanism and migrating on the membrane in different ways (see [10, 4, 34]). It is natural 
to ask if the Fleming- Viot convergence is valid in this general framework. The results in this 

^This convergence was proved in [19] using the technique of particle representation described in [!)]. This 
technique cannot be easily extended to the multi-type setting of this paper. Therefore the convergence proof 
in this paper is vastly different. 
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paper show that it is indeed the case as long as certain basic elements of the dynamics are 
preserved. This ensures that the analysis in [19] can be extended to more complicated (and 
realistic !) models for cell polarity. We discuss this example further in Section 3.2. 

Note that the geographical space E can be considered as the space of genetic traits. 
This casts our models into the setting of population genetics. The spatial migration can be 
seen as mutation that may happen at any time during the life of an individual, while the 
offspring dispersal mechanism is like mutation that can only happen at the time of birth of 
an offspring. We assume that the reproduction is clonal in the absence of mutation. The 
position-dependent birth and death mechanism is analogous to the selection mechanism in 
population genetics. Hence it is not surprising that spatial migration, offspring dispersal and 
position-dependent birth and death mechanisms correspond to mutation and selection in the 
limiting Fleming- Viot process. What is more interesting is that the sampling mechanism 
arises naturally from our models in the infinite population limit. This sampling mechanism 
is a key feature of the standard models in population genetics, such as the Wright-Fisher 
model, the Moran model and their variants (see [40, 29, 14]). This mechanism makes the 
models tractable by keeping the population size constant. It is done by matching the birth 
of an individual with the death of another individual chosen uniformly from the population. 
It is obvious that such a mechanism is quite unrealistic, at least for finite populations which 
are naturally fluctuating. However our Fleming- Viot convergence result shows that one 
can recover this sampling mechanism in the infinite population limit if the dynamics has 
an equilibrating density regulation mechanism that acts at a faster timescale than other 
events. It is well-known that a Fleming- Viot process arises in the infinite population limit 
of an appropriately scaled version of the Wright-Fisher or the Moran model (see [17] and 
[13]). Therefore if all the individuals have the same type (that is, q = 1) and heq is the 
equilibrium population density, then for a large (but finite) value of the scaling parameter 
A^, our models will have roughly the same dynamical behaviour as a suitably chosen Wright- 
Fisher or Moran model with the constant population size Nheq. This insight provides a 
justification for the assumption of a constant population size in population genetics models. 
Most of the mathematical literature on population genetics is concerned with two types of 
questions. In the absence of mutation, one wants to know the probability and the time of 
fixation of a particular genetic trait. The term fixation describes the event in which the 
whole population has the same genetic trait. In the presence of mutation, one attempts 
to investigate the properties of the stationary distribution, if such a distribution is present. 
These questions are difficult to answer for finite populations and one typically answers them 
by studying the limiting Fleming- Viot process. Our discussion shows that for large N, the 
fixation times and probabilities or the stationary distribution will be approximately the same 
for our model and the corresponding Wright-Fisher or Moran model. We illustrate this point 
through an example in Section 3.1. 

This paper is organized as follows. In Section 2 we describe the mathematical models 
that we consider and state our main results. In Section 3 we discuss the aforementioned 
applications of our results in greater detail. Finally in Section 4 we prove the main results. 
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Notation 



We now introduce some notation that we will use throughout this paper. Let M, M+, M*, 
N and No denote the sets of all reals, nonnegative reals, positive reals, positive integers and 
nonnegative integers respectively. For any a, 6 G M, their minimum is given by a A 6. 

Let ||-|| and (-, ■) denote the standard Euclidean norm and inner product in M" for any 
n G N. Moreover for any v = {vi, . . . ,Vn) G M", the norms \\v\\^ and \\v\\^ are defined 

ll^'lli = ^"=1 I'^il and \\v\\^ = maxi<j<„ The vectors of all zeros and all ones in 
W are denoted by 0„ and 1„ respectively. Let M(n, n) be the space of all n x n matrices 
with real entries. For any M G M(?T,,n), the entry at the z-th row and the j-th column 
is indicated by Mij. Its infinity norm is defined as ||M||^ = maxi<j<„ ^"^-^ \ Mij\ and its 
transpose and inverse are indicated by M'^ and respectively. The symbol /„ refers to 
the identity matrix in M(n, n). For any v = {vi, . . . , Vn) G M", Diag(w) refers to the matrix 
in M.{n,n) whose non-diagonal entries are all and whose diagonal entries are vi,. . . ,Vn- 
A matrix in M(?t,, n) is called stable if all its eigenvalues have strictly negative real parts. 
While multiplying a matrix with a vector we always regard the vector as a column vector. 

Let [/ C M" and V C M'". Then for any k G No, the class C^{U,V) refers to the set 
of all those functions / that are defined on some open set O C M" containing U such that 
f{x) & V for all X G t/ and / is fc-times continuously differentiable at any x G O. 

Let {S,d) be a metric space. Then by B{S) {C{S)) we refer to the set of all bounded 
(continuous) real- valued Borel measurable functions. If S is compact, then C{S) C B{S) 
and both B{S) and C{S) are Banach spaces under the sup norm \\f\\^ = sup^-g^ 
Recall that a class of functions in B{S) is called an algebra if it is closed under finite sums 
and products. Let B{S) be the Borel sigma field on S. By AipiS) and V{S) we denote the 
space of all finite positive Borel measures and the space of all Borel probability measures 
respectively. These measure spaces are equipped with the weak topology. For any / G B{S) 
and jj, G MpiS) let 



The space of cadlag functions (that is, right continuous functions with left limits) from 
[0, oo) to S is denoted by Ds[0,oo) and it is endowed with the Skorohod topology (for 
details see Chapter 3, Ethier and Kurtz [12]). The space of continuous functions from [0, oo) 
to S is denoted by Cs[0,oo) and it is endowed with the topology of uniform convergence 
over compact sets. An operator A on B{S) is a linear mapping that maps any function 
in its domain T>[A) C B[S) to a function in B{S). The notion of the martingale problem 
associated to an operator A is introduced and developed in Chapter 4, Ethier and Kurtz 
[12]. In this paper, by a solution of the martingale problem for A we mean a measurable 
stochastic process X with paths in Ds[0, oo) such that for any / G 'C(^), 



is a martingale with respect to the filtration generated by X. For a given initial distribution 
IT G V{S), a solution X of the martingale problem for A is a solution of the martingale 
problem for {A,7i) if vr = PX(0)^^. If such a solution X exists uniquely for all n G V^S), 
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then we say that the martingale problem for A is well-posed. Additionally, we say that A is 
the generator of the process X. 

Throughout the paper =^ denotes convergence in distribution. 

2 Model descriptions and the main result 

Our first task is to describe the models that we consider in the paper. As mentioned in 
Section 1, we model a population which resides in some compact metric space E and in 
which the individuals have one of q possible types. We denote these types by elements in 
the set Q = {l,2,...,g}. We identify each individual located at x & E with the Dirac 
measure 6x, concentrated at x. Moreover each individual is assigned a mass of 1/N where 
iV G N is our scaling parameter. For any i & Q, the population of type i individuals can be 
represented by an atomic measure of the form 



where rij is the total number of type i individuals and x\, . . . , G E are their locations. 
Define the space of atomic measures scaled by as 



Note that J^iy a{E) C Aip^E), where Ai f{E) is the space of finite positive measures. Let 
Ai'j^ ^{E) and J^p{E) be the spaces formed by taking products of q copies of Aijy^aiE) and 
JVlp[E) respectively. Since for each i E Q, the type i population can be represented by a 
measure /ij G J^N,a{,E) , the entire population can be represented by a g-tuple of measures 



Let 1e denote the constant function in C{E) which maps each point in to 1. Define 
the density map H : A4.'^p{E) — )■ as the continuous function given by 

H{^i, ...,iJ,g) = {{1e, yUi), . . . , {Ie, /ig)) for any /i = (/ii, . . . , /ig) G M^piE). (2.2) 

We will refer to h = H{fi) as the density vector corresponding to yU G Ai''p{E). Note that if 
the population is represented by a /i G M.%^{E) and ii h = {hi, . . . , hg) is the corresponding 
density vector, then hi is just the total number of type i individuals divided by A^. The 
density vector h contains no information about the distribution of individuals on E. 

2.1 The type-dependent migration mechanism 

In our models, each individual of type i E Q will migrate according to an independent E- 
valued Markov process with generator Bi. We will assume that each operator Bi generates a 
Feller semigroup on C{E) (see Chapter 4 in Ethier and Kurtz [12]). Furthermore we assume 





(2.1) 



/i = (/ii,...,^g) eM%^, 
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that there is an algebra of functions Vq C C{E) which is dense in C{E), contains 1e and 
satisfies 

Vo C V{Bi) for all i E Q. (2.3) 

The martingale problem corresponding to each Bi is well-posed and any solution is a strong 
Markov process with sample paths in De[0,oo) (see Theorem 4.2.7 and Corollary 4.2.8 in 
[12]). 

We now formally describe how this type-dependent migration of individuals translates 
into the evolution of our population in the space Ai'j^^{E). For each n G N, define a space 
of atomic probability measures as 



k j=i 



and a class of continuous real-valued functions over V{E) by 

Co = I^Fiu) = f[{fi, v) : h,...,fmeVo and m G n| . (2.4) 

Suppose that u = (1/n) Yl]=i G 'Pn,a and F^u) = YYiLiifh ^) ^ Cq. For positive integers 
k < m, let be the set of onto functions from {1, . . . , m} to {1, ... ,k} and for any p G P™ 
and / = 1, . . . , let 

f!'\^)= n /^■(^)- (2.5) 

iep-i(0 

Then we can write 

m / n \ n m 

1=1 \j=i / ii,...,j™=i 1=1 

m n k 

= -"'"EE E Wi"'^^nl (2.6) 
k=i p£p;j' ju...,jk=i 1=1 

where the last term has summation over distinct choices of ji, . . . , jk G {1, . . . , n}. For each 
i E Q,n eN we now define an operator B" : P(B[^) = Cq ^ B (Vn,a) by 

m n k k 

B-F{v)=n-Y.Y. E Y.B^fi'^^'^n) H fr'\^^rl (2.7) 

k=l peP^ ii,...jfc=l (=1 r=l,rj^l 

where F G Co is given by (2.6). Observe that any F G Cq is bounded and 

sup sup |B^F(z/)| < oo. (2.8) 

ngN ueVn^a 
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One can easily verify that the martingale problem for each B" is well-posed. If uq = 
(1/n) Yl]=i^xj € 'Pn,a then the solution of the martingale problem for (B",5i,(,) is just the 
empirical measure process of a system of n individuals moving in E according to independent 
Markov processes with generators Bi and initial positions xi, . . . ,Xn- For more details see 
Section 2.2 in Dawson[7]. 

For any /i, ...,/„, G I^o consider a function F : A^^(£') -> M of the form 

m / q \ 
F{f^) = ll[j2c,{h){f„^,)\, (2.9) 
j=i \i=i J 

where h = H{fi) and c : Mj^ — )■ is a function that satisfies 

sup (c(/i),/i) < oo. (2.10) 



Such a function F is bounded because 



ti&MliE) \l=l,...,m 



sup |F(/i)| < max ||/i||oo sup (c(/i), /i) 



+ 



Define a class of functions by 

= |f G 5 {M%{E)) : F is given by (2.9) and c satisfies (2.10)} . (2.11) 

If /i = (/ii, . . . , /iq) G A^^^(-E), then for each i E Q, fii has the form /ij = (1/A^) X]j=i ^x'- 
Let z/j = (l/?^i) X]j=i "^i:* if rij > and i^j = 6xo if = 0, where xq E E is arbitrary. Let 
be as before. Pick a F G Cg of the form (2.9). For any p G and / = 1, . . . , define 



-F/^"* G Co by fI^\v) = {fi^\ ly), where //^■' is given by (2.5). We can write the function F 

q m q m 

n,...,im = l i = l il,...,jm = l i = i 

fc=ipePfc"'ii,...A=i r-=i 

where the last term has summation over distinct choices of G {l,...,g}. Let 

gAT . 2)(B^) = 5 {M%SE)) be the operator whose action on any F G Cq written in 

this form is given by 



as 



B^F(/x) = E E E (Q.(^)^u)'^"^^^' ■ ■ ■ Mh)hi 
k=ip&P]^h,...,ik=i 

k k 

x^B^F(-)(^J n ^^'VU, (2.12) 
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where for any n G N, i G the operator B" is given by (2.7). For convenience is defined 
to be the identity map on Cq. The function B^F is bounded due to (2.8) and (2.10). 

The martingale problem for B^ is well-posed because the martingale problem for B" is 
well-posed for each i G Q,n G N. To see this suppose that fio = (/io,i, • • • , /^o,.?) ^ -^n ai^) 
and for each i E Q, /io,i G M.j~^a{E) has the form /io,j = N)YTj=i^x^-- If 'T-j > let 
z/Qi = iX/n-i) YTjLi ^x^. ci-iid if rij = let i^o,i = ^xo for some arbitrary xq G E. For each i E Q, 
let {z/j(t) : t > 0} be the unique solution of the martingale problem for (B"%(5yg .). Then the 
process {fi(t) : t > 0} given by 

Mt)=(^^i(t),...,^^,(t)), 

is the unique solution to the martingale problem for (B^,(5^q). 



2.2 The density regulation mechanism 

We mentioned in Section 1 that in our models, the birth and death rates of individuals 
depend on the population density in such a way that they induce an equilibrating density 
regulation mechanism. We now describe this mechanism in our multi-type setting. 

For each i,j G Q, let Pij be a bounded function in C^(R^_,M_|_) and for each z G Q let 
Pi be any function in C^(]R^,M+). For any h G M^, let the matrix A{h) G M(g, g) and the 
vector 9{h) G M'^ be given by 

Mh) = i , ..t^^^M r""' 1' (2-13) 
^ \ Pii{h) - pi{h) otherwise J' ^ ^ 



and 



9{h) = A{h)h. (2.14) 



Consider 6' to be a vector field over M^. Observe that if G is such that hi = then 
6i{h) > 0. This shows that any solution to the initial value problem 

^ = e{h), hiO) = hoeRl (2.15) 

stays inside for all positive times for which it is defined. Standard existence and unique- 
ness theorems imply that for any Hq G M!j_, there is a solution h{t) of (2.15) defined on some 
maximal time interval [0,a). Moreover if a < oo then ||/i(t)||]^ — ?■ oo as t — ?■ a_. Since is 
bounded for each i,j G Q, there is a positive constant Cq such that 

1 Q 

Xl^i(^) ^ J2 f^i'W^j ^ Ce\\h\\i for all /i G (2.16) 

i=l i,j=l 



and hence 



d\\hm, _^dh,{t) 



J2!^ = Y,0^^h{t))<Ce\\hm,. (2.17) 



dt ^ dt 

«=1 2=1 
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Therefore using Gronwall's inequality and (2.17) we obtain that < ||/i(0)||j^ e*"** for 

all t G [0, a) and so cannot go to oo as t — )■ a_. Thus a = oo and this shows that for 

any ho G M+, the initial value problem (2.15) has a unique solution which is defined for all 
t > 0. 

Let ipe '■ X M+ — Mj^ be the flow associated to the vector field 9. This means that ipe 
satisfies 

ipe{x,t) = x+ [ e{'4}g{x,s))ds for all (x,t) G X R+. (2.18) 
Jo 

This fiow is well-defined because of the arguments given in the preceding paragraph. In fact 
since 9 is in C2(]R9,M3), the map tjje is in C2(]R« x ]R+,M!^). This map also satisfies the 
semigroup property 

i/je{x,t + s) = i/jg{ilJe{x,t),s) for all x G M!^ and s, t G M+. (2.19) 

We will say that a set f/ C M+ is ^/'e-invariant if for all t > 0, ipg{U, t) G U where ipeiU, t) = 
{ipg{x,t) : X G f/}. Before we proceed, we need to make some more assumptions. 

Assumption 2.1 (A) There exists a vector h^q G M+, h^q ^ 0^ such that 9{heq) = 

A(h(,q)hcq = Oq. 

(B) The Jacohian matrix [J9{hcq)\ is stable. 

(C) The matrix A{hcq) is irreducible, that is, there does not exist a permutation matrix 
P G M(g, q) such that the matrix PA{hcq)P~^ is block upper-triangular. 

(D) For each i,j G Q, the functions pi and are analytic at h^q, that is, they agree with 
their Taylor series expansion in a neighbourhood of heq. 

Part (A) says that there is a nonzero vector hcq G which is a fixed point for the fiow 
ipg. Part (B) implies that this fixed point hcq is asymptotically stable for this fiow. The 
significance of part (C) will become clear later in this section. Part (D) is a technical 
condition that we require to prove our main result. 

We define the region of attraction of the fixed point h^q for the fiow ipg by 

Ucq=\heRl: lim Mh, t) = hA • (2.20) 

I. t— s-oo J 

Part (B) of Assumption 2.1 and Lemma 3.2 in Khalil [21], ensure that Ucq is a '0£)-invariant 
open set in Note that Ucq may not be an open set in M'^. 

We are now ready to describe the density regulation mechanism. Suppose that when 
the population density vector is /i G M+, then for each i,j G Q, an individual of type i 
gives birth to an individual of type j at rate Pijih) and an individual of type i dies at rate 
Pi{h). At the time of birth, the offspring is placed at the same location in E as its parent. 
Note that the birth and death rates do not depend on the location of the individuals. If the 
scaling parameter is G N and the mass of each individual is 1/A, then we can view this 
density- dependent population dynamics as a Markov process over the state space A^^^(-E) 
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with generator : V(K^) = B{M%^{E)) ^ B{M% ^{E)) defined as follows. For any 
F e B{M%SE)) and any /i e -M^,„(^) with h = H{fj.) ' 

R^F(/.) = J2 Njj.,{h) (/i + j^6i^ - (2.21) 



[ p^{h) (f L - Ui) - Fif,)) 



where for any = (/xi, . . . , /i^) G A^^^(-E'), j E Q and x E E 

Concrete results on the well-posedness of the martingale problem for R'^ will come later. 
First let us understand how this operator R^ drives the population density to an equilibrium 
value. Let {j^^it) : t > 0} be a Al^^(E)-valued Markov process with generator R^ and 
let {h^(t) : t > 0} be the corresponding density process defined by h^(t) = H{fi'^{t)). 
Assume that h^{0) — )■ ho a.s. as — oo and is some vector in Ucq- From Theorem 
11.2.1 in [12], one can conclude that a.s N oo, {h^{t) : t > 0} converges in the Skorohod 
topology in D]Rq[0,oo) to the process {ipeiliQ.t) : t > 0}. Since G f/cq, ipeih^.t) — > /icq as 
t — )■ oo which indicates that for a large A^, the density process /i^(-) gets closer and closer to 
/icq with time. This shows, at least informally, that the operator R^ drives the population 
density towards h^^. Henceforth we shall refer to the vector /igq as the equilibrium population 
density. Of course, this discussion totally ignores how R^ affects the spatial configuration 
of the population. Our main goal in this paper is to discover how the spatial distribution of 
the population evolves when the density is equilibrated at a faster timescale in comparison 
to the other mechanisms. 

For any h G M+, the matrix A{h) given by (2.13) signifies how the various types of 
individuals interact when the population density is h. Part (C) of Assumption 2.1 means 
that at the equilibrium population density, all the types of individuals are communicating 
with each other by influencing each other's birth and death rates. 



2.3 Mathematical Models 

We now describe the models that we consider in this paper. We start with a basic model 
which only has spatial migration along with density regulation. We then extend this model 
by adding other features such as position dependence in the birth and death rates, offspring 
dispersal and immigration. All the models will be parameterized by the scaling parameter 
A^, with 1/A^ being the mass of each individual in the population. We will describe each 
model by specifying the generator of the associated Markov process. The well-posedness of 
the martingale problems corresponding to these generators is given by Proposition 2.4. Our 
main results are presented in Section 2.4. 
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2.3.1 Basic Model 

In this model, the individuals migrate according to the type-dependent migration mechanism 
specified in Section 2.1 and their birth and death rates regulate the population density as 
described in Section 2.2. If the scaling parameter is A^, then at any time we represent the 
population as a measure in A^^^(£'). The population evolves due to the following events. 

• Each individual of type i E Q migrates in E according to an independent Markov 
process with generator Bi. 

• When the population density vector is /i G M+, each individual of type i E Q gives 
birth to an individual of type j E Q at rate NPij[h). At the time of birth, the offspring 
is placed at the same location as its parent. 

• When the population density vector is h E M+, each individual of type i E Q dies at 
rate N pi{h). 

This population dynamics can be viewed as a A^^^(-E')- valued Markov process whose gen- 
erator : V{K^) =Cl^ B{M%^^{E)) is given by 




(2.22) 



for any F E Cq. Here h = H{fi) is the density vector corresponding to /i. 



2.3.2 Model with position dependence in the birth and death rates 

In the above model, the birth and death rates of individuals do not depend on their loca- 
tion. One may want to consider models in ecology where some spatial locations are more 
advantageous for reproduction or some locations are more hazardous for survival. If we think 
of E as the space of genetic traits, then one may consider models in which the trait of an 
individual infiuences its chances of reproduction or survival. To capture such situations we 
now introduce a model in which the birth and death rates of an individual can vary with 
its position. However we will assume that this position dependent variation is small, in the 
sense that even though an individual's birth and death rate is of order (as in Section 
2.3.1), the spatial variation in these rates is of order 1. We now make the model precise. 

For each i,j E Q, let b^j, be bounded continuous functions from E x M5j_ to M+. These 
functions determine the spatial variation in the birth and death rates. The migration of 
individuals is like in the basic model. However now the birth and death mechanism changes 
as follows. 

• When the population density vector is h E M^, an individual of type i E Q located at 
X E E, gives birth to an individual of type j E Q at rate bfj{x, h) + N[3ij{h). At the 
time of birth, the offspring is placed at the same location as its parent. 
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• When the population density vector is /i G M^, an individual of type i E Q located at 
X E E, dies at rate h) + Npi{h). 

The evolution of our population under this dynamics can be viewed as a A^^^(i?)-valued 
Markov process with generator : V{Af) = Cq — )■ i?(A^^^(£')) defined for any F G Cg 
by 



A^Fif^) = B''F{f,) + J2 N J^{bUx,h) + N(3,,{h)) (^F (/^ + ^5^) - ^^^{dx) 

(2.23) 

+ [ idtix,h) + Np,ih)) (f L-hl]-Fip)]fi,idx), 

i&Q V V / / 

where h = H{fi) is the density vector corresponding to fi. 
2.3.3 Model with offspring dispersal at birth 

In the basic model we described in Section 2.3.1, the offspring is placed at the same location 
as its parent at the time of birth. However we may want to construct models where this 
restriction needs to be relaxed. For example, while modeling plant populations, one may 
wish to account for the spreading of seeds due to wind and other factors. Also in models 
for population genetics, where E is the space of genetic traits, offsprings may be born with 
a different trait than their parents due to mutations. To consider such situations we now 
present a model in which an offspring may be born away from its parent. We allow this 
offspring dispersal to either be rare (happens with probability proportional to 1/A^) or small 
(the offspring is placed at a distance proportional to 1/N from the parent). We handle both 
these cases in a unified way. 

For each i,j G Q and N E N, let {Ifj be a function from E to V{E) and let pfj be a 
function from E to [0, 1]. The individuals migrate and die in the same way as described in 
the basic model (Section 2.3.1). The birth rates are also the same as in the basic model. 
However when an individual of any type i G Q located at x G gives birth to an individual 
of type j, the location of the offspring is x with probability (1 — pfj{x)) and distributed 
according to i!^fj{x, ■) with probability pfj{x). 

To pass to the limit — oo, we need an assumption on pA and "dfj which is stated 
below. 

Assumption 2.2 For each i,j G Q we assume that there is an operator Cij : I){Cij) — )■ 
C{E), whose domain is taken to be the same as Vq (see (2.3) j for convenience, such that 
for every f eVq. 



lim sup 



Af-j-oo 



NpfAx) / {f{y)-f{x))dUx,dy)-C.,f{x) 



E 



0. 



The evolution of our population under the dynamics described above can be viewed as a 
A^^_„(E)-valued Markov process with generator A^ : ^'(A^) = Cl ^ B{M%^^{E)) defined 
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by its action on any F G Cg by 

A^Fi^^) = B^F(/.) + N'J2 [ P^W (f L - hi) - (2.24) 



fii{dx), 



where h = H{fi) is the density vector corresponding to fi. 
2.3.4 Model with immigration 

Consider a population whose dynamics is as described in the basic model (Section 2.3.1). In 
addition, suppose that the individuals of each type are immigrating to at a certain density 
dependent rate and settling down according to some distribution on E. In this section we 
model this situation. Such a model can help us understand the effects of immigration on the 
population demography. 



For each i E Q,\et Ki 



be a continuous function satisfying 



^) for all h G 



(2.25) 



for some C > 0. The individuals migrate, reproduce and die as in the basic model. Moreover, 
when the population density vector is /i G M+, the individuals of each type i E Q arrive in the 
population at rate Nni{h) and their initial location is given by the distribution Gj G 'P{E). 

The evolution of our population under this dynamics can be viewed as a A^^^(-E')-valued 
Markov process with generator : V^A^) = Cg — ?■ B{A4.'^^ ^{E)) defined by its action on 
any F G Cg as 



(2.26) 



F(^) = B^F(/i) + «:,(/z) / [fL+ - F{f,)] 0,(dx) 

i,j^Q ^ ^ ^ ^ 

where h = H{fi) is the density vector corresponding to /x. 

Remark 2.3 For each / G {0,1,2,3} define the operators Gf : V{Gf) = B{M%{E)) -> 
B^M^p^E)) as follows. For any F G 5 {M%{E)) let 



G^F{^^) =0, 



Gf F(/i) =N 



J2 1^ h) (^F (^/i + ^6i}j - F(/i)^ ^,(dx) 



(2.27) 
(2.28) 
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and GgM/i) =iVE'^^(^) / U + l^^i) - Hf^)) 0' 

iGQ Je \ \ J J 



[dx) 



where h = H(fi). Then for each / G {0, 1, 2, 3} and F E Cq we can write 

Af F = B^F + iVR^F + F. 



(2.29) 
(2.30) 

(2.31) 



This form makes the timescale separation clear between the fast density regulation mech- 
anism (A^R^) and the slow migration (B^), position-dependent birth and death (Gf'), 
offspring dispersal (G^) and immigration (G^) mechanisms. □ 



2.4 The main results 

The main results of this paper are concerned with the limiting behaviour of the dynamics 
under the models described in Section 2.3. Before we present these results we must first 
verify that all the models in Section 2.3 can be represented by a suitable Markov process. 
This is established by the following proposition which will be proved in Section 4.1. 

Proposition 2.4 For each I G {0, 1, 2, 3} and N E N, the F>j^<i^ {£;)[0, oo) martingale prob- 
lem for Af is well-posed. 

We now begin analyzing how a sequence of Markov processes with generators Af^ behave 
as iV — oo. The next proposition exhibits some important properties about the limiting 
dynamics. The proof of this proposition is given in Section 4.2. Recall the definition of Ucq 
from (2.20). 

Proposition 2.5 Fix a / G {0, 1,2,3}. For each N e N, let {/i^(t) : t > 0} be a solution 
of the martingale problem for Af and let {h^{t) = H{^^{t)) : t > 0} be the corresponding 
density process. Assume that there is a compact set Kq C Ucq such that h^{0) G Kq a.s. for 
all G N. Let t^v be a sequence of positive numbers satisfying — > and Nt^- oo as 
N ^ oo. Then we have: 

(A ) For allT >0 

sup Wh'^ {t + tN) — hcq\\-^^ ^ as N ^ oo. 
te[o,T] 

(B) For allT>0, f E C{E) and i,j eQ 

sup \hf{t + tN){f,fi^{t + tN))-hf{t + tN){f,fif{t + tN)) \ ^0 asN^oo. 
te[o,r] 
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Let the processes : t > 0} and : t > 0} be as in the above proposition. 

Part (A) of this proposition imphes that the process h'^{- + tjy) =^ h^q as — > oo in 
Dk<;[0, oo). In other words, for large A^, the density process is constantly near the equilibrium 
population density /leq (after a small time shift In). This emphasizes the point that we 
made in Section 1. The density regulation mechanism operating at a faster timescale than 
our timescale of observation, keeps the population density equilibrated at all times. Note 
that the process {/i^(t) = (/i^(t), . . . , fi^ (t)) : t > 0} is A/l^(-E) -valued and it keeps track 
of how the populations corresponding to all the g-types are evolving in the space E. Part 
(B) of the above proposition shows that in the limit, all the q sub-populations are spatially 
fused (in proportions determined by the density vector). Hence their spatial evolution can be 
studied together by using a single 'P(£')-valued process. This kind of model reduction result 
is quite common in stochastic reaction networks with multiple timescales (see [3] and [22]), 
where one can often equilibrate the concentrations of the fast chemical species and derive a 
reduced model for the dynamics of the slow species. In our case, the dynamics at the fast 
timescale equilibrates the population density as well as the relative abundances of all the q 
sub-populations at each location on E. These relative abundances equilibrate because the 
birth-death interaction matrix A{h) (see (2.13)) is irreducible at the equilibrium density heq 
(see part (C) of Assumption 2.1). 

The proof of Proposition 2.5 will exploit the form (2.31) of the operator Af^ . A brief 
outline of the proof is as follows. We will define a M.\. x valued process {X^{t) : t > 0} 
by 

X^(t) = (hnt), . . . , C(t), Y.'^it), . . . , Yf_,{t)), (2.32) 

where each Y-^{t) is a density-dependent linear combination of terms like {h^ (t) {f , fif (t)) 
—hf{t){f, it))) for some choice of / G C{E). Next we will show that is a semimartin- 
gale which satisfies an equation of the form 

X^(t) = A:^(0) + N f F(X^(s))rfs + Z^(t), (2.33) 

Jo 

where {Z^ : A^ G N} is a sequence of ]R^'^~^-semimartingales which is tight in the space 
D^2q-i[0, oo). This clearly indicates that for large values of A^, the drift term of the form 
NF{X^ {■)) completely overwhelms the effect of the semimartingale Z^. Equations like 
(2.33) were studied by Katzenberger in [23] in a much more general setting. He showed 
that under certain conditions, the sequence of semimartingales {X^ : A^ G N} converges in 
distribution to a semimartingale X as A^ — )■ oo. Moreover X only takes values in a set F 
which is an invariant manifold for the deterministic flow induced by F. In our case, this set 
r only consists of one point Xcq = (/icq, Oq~i) and this enables us to prove Proposition 2.5. 
The details are given in Section 4.2. 

We mentioned before that in the limit A^ — )■ oo, the spatial evolution of all the q sub- 
populations is governed by a single V{E)-valued process. Our next result. Theorem 2.6, 
shows that this ■p(£')-valued process is in fact a Fleming- Viot process that can be character- 
ized by its generator. Before we state Theorem 2.6 we first need to introduce several objects. 
The existence and properties of some of these objects will be studied in the appendix. 

Recall the equilibrium population density vector heq = (/ieq,i, • • • , ^eq,(j) from Section 
2.2. It can be verified that this vector has strictly positive components (see part (A) of 
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Lemma A.l). Moreover part (C) of Lemma A.l shows that there is a unique vector Veq 

(Weq,l, . . . , Ueq,g) ^ SUCh that 



A{hcq) = Oq and (Wcq, /icq) = ^Veq,iheq,i = 1- (2.34) 



1=1 



Observe that Vq C C{E) satisfies (2.3). Define an operator i?avg : ^^(-Bavg) =1^0—^ C{E) by 

Bi,^gf = ^Vcq,ihcq,iBif for f eVo. (2.35) 

From (2.34) we can see that the operator i?avg is a convex combination of the operators 
{Bi : i G Q}. Let 7smpi be the positive constant given by 

7smpl = vlq,Aq,iPiihcq)- (2.36) 

For each i,j G Q, let the functions fof^, (i| be as in Section 2.3.2. Define fe^vgi '^avg ^ C'(-E') as 

Kvgi^) = Y bij{x,heq)Veq,jheq,i and <vg(x) = ^ (i- (x, /leq)t'cq,i/loq,i for X G -E. 

For each z,j G Q, let the operator Cij be as in Assumption 2.2. Define the operator Cavg : 
^^(Cavg) =Vo^ C{E) by 

Ca.vgf = (3ij{hcq)Vcqjhcq^ ,a,f for / G Po- (2.37) 

For each i E Q, let Kj,9j be as in Section 2.3.4. Define the operator Javg : 'D(/avg) = 
B{V{E)) ^ B{V{E)) by 

/avg/(x) = ^ ^^^{K^Kq, ! (/(y) - /(x))0,(rfy) for / G i?(P(i?)). (2.38) 

We now define the operators Aq, Ai,A2 and A3 with domain Co (see (2.4)) as below. For 
anyF(i.) = n;^,(/,,z/)GColet 

m 

A,F{v) = Y,{B..Ji, 1^) U^fv ^) (2.39) 

~l~ O'smpl 

l<l=ik<m j=a,k 
m 

1=1 j^i 

m. 

1=1 



m 



(2.41) 



1=1 j^i 



m 



and AsFiu) = A,F{u) + ^(/avg/«, y) \{{f,. v)- 



(2.42) 



We will assume that the operators -Bavg, (-Bavg + C'avg) and (-Bavg + -^avg) generate Feller semi- 
groups on C{E). The well-posedness of the martingale problems corresponding to Aq, Ai,A2 
and A3 follows from Theorem 3.2 in [13]. In fact, any solution will have sample paths in 
C-p(s) [0, 00). The operator Aq is the generator of a neutral Fleming- Viot process on E 
with mutation operator i?avg and sampling rate 27smpi- The operators A2 and A3 generate 
a similar Fleming- Viot process with the mutation operator changed to (i?avg + Cavg) and 
(i?avg + -^avg) respectively. The operator Ai also generates a similar Fleming- Viot process, 
but with selection. The last two terms in its definition correspond to fecundity selection 
(with intensity function bl^^) and viability selection (with intensity function rf^^g). See Don- 
nelly and Kurtz [8] for more details. We now formally state the main result of our paper. 
The proof is given in Section 4.4. 

Theorem 2.6 Fix a I G {0,1,2,3} and let {fi^{t) : t > 0} be a solution to the martingale 
problem for A^. Suppose that /i^(0) =^ /i(0) as — t- 00 and i7(/x(0)) G f/cq a-s. where 
Ucq is given by (2.20). Let t^- be a sequence as in Proposition 2.5. Define another process 



Then there exists a distribution ir G V {V{E)) such that =^ heqi' in Z}_vi|,(£;)[0, 00) as 
iV — >■ 00 and {I'it) : t > 0} is a Fleming-Viot process with type space E, generator A/ and 
initial distribution n. 

Remark 2.7 The initial distribution n of the process {i^(t) : t > 0} is related to the distri- 
bution of fi{0). This relation is stated in Remark 4-10. 

Remark 2.8 In Section 2.3 we first defined a basic model and then constructed auxiliary 
models by adding other mechanisms, one at a time. These other mechanisms are position 
dependent birth and death, offspring dispersal and immigration. One can consider models 
in which more than one of these mechanisms are simultaneously added to the basic model. 
The proof will demonstrate that the generator of the limiting Fleming-Viot process is then 
obtained by adding the correct term corresponding to each of these additional mechanisms to 
the operator Aq. This correct term can be seen from the definitions of Ai, A2 and A3. For 
example, one can have the basic model along with position dependent birth and death (Section 
2.3.2) and offspring dispersal (Section 2.3.3). Then the limiting Fleming-Viot process has 
the generator given by 



{/2^(t) ■.t>0}by 



^^(t) =^^(t + tjv) forty 0. 



(2.43) 



m 
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i=\ jj^i 1=1 j^i 

for any F{u) = UZiifi^ e Cq. 

We now give a heuristic explanation of why the dynamics under the models described in 
Section 2.3 converges to a Fleming- Viot process. Note that part (A) of Proposition 2.5 says 
that the population density is pinned to a constant value h^q in the limit. Therefore any 
addition of new mass in the population must be concurrently offset by an equal reduction 
of existing mass and vice versa. Furthermore, when the mass is reduced or added to keep 
the balance, this reduction or addition happens at locations that are chosen more or less 
uniformly from the current empirical measure of the population. This is because the birth 
and deaths rates of individuals are dominated by a term which is density dependent but 
location independent. This argument offers some intuition as to why the fast birth-death 
terms (that form part of the operator A^R^) give rise to the sampling term in the limit (the 
second term in Aq). It also shows why the position dependent birth and death terms in A^ 
become selection terms in Ai and the offspring dispersal (immigration) term in A^ (A^) 
becomes a mutation term in A2 (A3). Since the position of an individual in E can also be 
seen as its genetic trait, one can interpret the migration on E as genetic mutation. Hence it 
is not surprising that the migration operators appear as part of the mutation operator in the 
limiting process. Part (B) of Proposition 2.5 says that in the limit, all the q sub-populations 
become spatially inseparable. This causes all the mechanisms in the limiting process to 
appear in an averaged form. 

Let {/i^(t) = {fii{t), . . . , fig (t)) : t > 0} be a Markov process with generator Af, for 
some I e {0, 1, 2, 3}. It is difficult to prove the convergence of this process directly because 
the density regulation mechanism acts on it at the fast timescale. This can be seen by 
splitting the operator A^ according to (2.31) and noting that A^R^ becomes unbounded as 
N 00. To pass to the limit we consider another measure- valued process {z/^(t) : t > 0} 
that is constructed by suitably combining the various components of {fi^{t) : t > 0}. In 
particular 

z/^(t) = ^ A.(/i^(t))/if (t) for t > 0, (2.44) 

i=l 

where h^{t) = H{fi^{t)) and A = (Ai, . . . , A^) is a function from to which satisfies 
certain conditions. These conditions are chosen to ensure that {z/^(t) : t > 0} is a V{E)- 
valued process whose dynamics is such that the density regulation mechanism acts at the 
slow timescale. Such a function A can be shown to exist by proving that a certain system 
of coupled partial differential equations has a solution with some desired properties. This 
is done in Section 4.3. We will then show that as A^ — )■ 00 we have ^ v where {^{t) : 
t > 0} is the Fleming- Viot process specified by Theorem 2.6. This convergence along with 
Proposition 2.5 allow us to prove Theorem 2.6. The details of the proof are given in Section 
4.4. 

The discussion in the preceding paragraph also shows that intuitively we can think of 
the limiting Fleming- Viot process as describing the spatial evolution of a mixed popula- 
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tion formed by taking a suitable density-dependent linear combination of all the q sub- 
populations. This is reminiscent of the notion of virtual species (formed by linearly combin- 
ing several chemical species), that are needed in the specification of the reduced models in 
chemical reaction networks with multiple timescales (see [5]). 



3 Applications 

In this section we discuss the applications mentioned in Section 1 in greater detail. Note 
that a Fleming- Viot process usually has continuous paths (see [13]). Hence Theorem 2.6 can 
be seen as a diffusion approximation result which shows that a stochastic process with jumps 
can be approximated by a process with continuous paths. Such results provide a justification 
for drawing inferences about the original process (with jumps) by analyzing a more tractable 
process with continuous paths. 

To demonstrate the usefulness of Theorem 2.6 we present two examples. In the first 
example we consider a population genetics model having logistic interactions along with 
rare mutation and weak selection. The words rare and weak indicate that the mutation 
and selection events occur at a slower timescale than other events. The difference between 
this model and a standard population genetics model (Wright-Fisher or Moran) is that 
the population size is not fixed but fluctuating due to the logistic interactions. Theorem 2.6 
guarantees that by taking the infinite population limit in a suitable way, we obtain a Fleming- 
Viot process. In many cases this limiting process is well-studied and using its properties one 
can estimate fixation probabilities, fixation times and the stationary distribution for the 
finite population model. Our second example sheds light on the phenomenon of cell polarity 
which refers to the spatial crowding of molecules on the cell membrane. We draw upon 
our work in [1*)] to show that Fleming- Viot convergence can help us understand how cells 
establish and maintain polarity. In [19] we only consider a very simple model, but the results 
in this paper ensure that the same analysis holds for a general class of models. 



3.1 Logistic model for population genetics 

The logistic growth model is very popular in ecology. It was proposed by Verhulst [37] in 
1838 to describe the growth of a population in the presence of competition for resources. In 
this model each individual reproduces at rate /3 and dies at a rate pP/N, where P is the 
current population size and is the carrying capacity of the habitat. In the deterministic 
setting, the population size (P) evolves as a function of time {t) according to the ordinary 
differential equation 

dP ^ P2 

Let h{t) = P{t)/N be the population density at time t. Then the above differential equation 
becomes 

^=f3h-ph^. (3.45) 
It is immediate that if h{0) > then h{t) — h^q := (3/p as t ^ oo. 
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We now construct a population genetics model that has logistic interactions along with 
rare mutation and weak selection. Suppose the compact metric space E is the set of all the 
genetic traits that an individual can have. Each individual is given a mass of with N 
being the carrying capacity as before. The population at time t can be represented by the 
measure 

/^''(^) = ]^ E (3-46) 

i=l 

where n^(t) is the number of individuals at time t and Xi,X2,--- G E are their genetic 
traits. Let b"^ be a continuous function from E to M+. When the population density (total 
mass) is h, an individual with trait x & E gives birth at rate {(3 + b^{x)/N) and dies 
at rate ph. Its offspring has the same trait x with probability (1 — p{x)/N). However 
with probability p{x)/N, the offspring is a mutant and its trait is chosen according to the 
distribution 'd{x,-) G V{E). The process {7I^(t) : t > 0} can be viewed as a Markov 
process with state space M.M,a{.E) (see (2.1)). The timescale at which we have described the 
dynamics is such that the mutation and selection events will vanish in the limit — )■ oo. 
Therefore to study their effects, we must observe the process at the timescale which is 
times slower. Let 

/i^(t) =ji^{Nt) fort>0. 

The dynamics of {/i^(t) : t > 0} has fast density regulation along with position-dependent 
birth (see Section 2.3.2) and offspring dispersal mechanism (see Section 2.3.3). Assuming 
that /i^(0) ^ ;u(0) as N — 7- OO and (l£;,yu(0)) > a.s. Theorem 2.6 gives us the following. 
If tjsi is a sequence satisfying — )■ and Nt]^ — oo, then the process /i^(- + tj^) ^ heqi^ as 
— )■ oo, where h^q = P / p and {^{t) : t > 0} is a Fleming- Viot process with generator given 

by 

AF{v) = PY.\I p(^) { I (/'(2/) - fi{xmx, dy)] u{dx)] \{{f,, v) (3.47) 

m 

1=1 j^i 

l<l^k<m j^l,k 

for any F{iy) = YYiLiifi^ ^) where fi, . . . , fm G B{E). This is of course the generator of a 
Fleming- Viot process with mutation and fecundity selection. We now present a couple of 
cases where this process is well-studied. 

Suppose that E = {!,... ,K} for some K E N. Then for any t > we can express 
z/(t) G V{E) as the K-tnple (i^i(t), . . . , i^i^(t)), where Uiit) is the proportion of individuals 
having genetic trait i E E. This representation allows us to view {iy{t) : t > 0} as a process 
over the i^-simplex 

r ^ 

= \ {xi, . . . , Xk) '■ Xi > and Xi = 1 
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For all i,i E E set 6ij = f3pii)'&{i, {j}) and = b'^{i). Then {iy(t) : t > 0} is a diffusion 
process over with generator given by 

where / G C^(]R^,M), u = {ui, . . . ,vk) and 5ij is the Kronecker delta function. This is 
the generator of the Wright-Fisher diffusion process [25]. Many explicit results about the 
fixation probabilities, fixation times and the stationary distribution can be found in [11]. 

Let us return to the situation where is a general compact metric space and the dynamics 
evolves according to (3.47). Assume that for all x G -E" we have h^{x) = 0, p{x) = 1 and 
'i9(x, ■) = "^oi'), for some non-atomic probability measure ■i^o ^ ^(-^)- The resulting Fleming- 
Viot process {z^(t) '■ t > 0} arises as a reformulation of the infinitely-many-neutral-alleles 
model due to Kimura and Crow [26] (see [12] and Section 9.2 in [13] for more details). In 
this case, i/(t) can be written as a countable sum ciiSxi for any t > (see Theorem 7.2 
in [13]). This means that at time t, ai fraction of the population is located at Xi. Arranging 
these aj-s in descending order we can extract a process over the ordered infinite simplex 

oo 



[xi, X2, . . . ) : Xi > X2 ■ ■ ■ > and Xi = 1 



i=l 



This extracted process is a diffusion process over Aoo whose various properties are presented 
in [11]. Furthermore in [13] it is shown that the Fleming- Viot process {I'it) : t > 0} is 
ergodic and its unique stationary distribution 11 G V(V{E)) is given by 

U{S) = P G for all S G B{V{E)), (3.48) 

where the infinite vector {(pi, (j)2, . . .) has the Poisson-Dirichlet distribution with parameter 
f3/2p and .^i, ^2, • • • are i.i.d. with distribution -do, independent of (0i, 02, . . . ). The Poisson- 
Dirichlet distribution was introduced and studied by Kingman [27] in 1975. 

The results mentioned in the last two paragraphs indicate the behaviour of the evolu- 
tionary dynamics under our original model for large values of N. 



3.2 Cell polarity 

Cell polarity is an important phenomenon and understanding the mechanisms responsible 
for it is a matter of fundamental concern for biologists. It is widely accepted that polarity 
is established in 3 stages (see [10, 2, 33]), which can be described as follows: 

1. An unpolarized cell receives a spatial cue that may be intrinsic (coming from inside 
the cell) or extrinsic (coming from the extracellular environment). 

2. This cue is interpreted by the membrane-bound receptor molecules. 
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3. The feedback network inside the cell is activated, which amplifies the weak initial signal 
into a robust signal that can direct the molecules towards the clustering site. 

The feedback network has two components : positive feedback which enables the membrane 
molecules to pull the cytosol molecules to their location on the membrane, and negative feed- 
back that pushes the membrane molecules into the cytosol. Positive feedback is responsible 
for the localized recruitment of molecules on the membrane while negative feedback helps 
in regulating the population size on the membrane. The molecules diffuse slowly on the 
membrane but rapidly in the cytosol. 

Even though the feedback mechanism may bring the molecules together on the membrane, 
any clusters that form may not persist due to spatial diffusion. This caused some biologists 
to propose that other additional mechanisms are needed to generate spatial asymmetry (see 
[18, 20]), but these mechanisms are not always found in cells that exhibit polarity. Hence 
it is important to investigate if the feedback mechanism can alone counter spatial diffusion 
to establish cell polarity. For this purpose, Altschuler et. al. [2] formulated a simple model 
based on the mechanisms mentioned above. We now describe their model. Consider the cell 
to be a sphere of radius i? in M'^. The whole cell has N molecules which may be present on 
the membrane or in the cytosol. The following four mechanisms change the configuration of 
molecules in the cell. 

1. Association mechanism: Each molecule in the cytosol can move to a uniformly chosen 
location on the membrane at rate kon- 

2. Positive feedback: Each molecule on the membrane pulls another molecule from the 
cytosol to its location at rate /cfb x (fraction of molecules in the cytosol). 

3. Negative feedback: Each molecule on the membrane is pushed into the cytosol at rate 

4. Spatial migration: Each membrane molecule is constantly diffusing on the membrane 
according to an independent Brownian motion with diffusion rate D. 

The association mechanism provides the initial spatial cue to trigger cluster formation. In 
[2], this spatial cue is intrinsic because the authors are concerned with spontaneous cell 
polarity, which means that polarity is established without any extracellular influence. Hence 
the association mechanism acts uniformly on the membrane. When one wants to consider 
polarity that is established in response to a chemical gradient (see [39]) then a molecule 
associating itself to the membrane must choose its location according to some distribution 
that encodes the gradient information. We mentioned in Section 1, that the positive feedback 
mechanism is like a birth process, where the pulled cytosol molecule is the offspring of 
the recruiting membrane molecule. This introduces genealogical relationships between the 
membrane molecules. A set of membrane molecules are said to belong to a clan if they have 
a common ancestor. Note that when the diffusion rate {D) is small, we would expect the 
clan members to be huddled together. 

The analysis of the above model in [2] gives some interesting results. When the dynamics 
is described deterministically, using a reaction-diffusion partial differential equation, then the 
model fails to capture cell polarity. However in the stochastic setting, the model does predict 
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the formation of clusters in certain parameter regimes, when the number of molecules (A^) 
is small. This result is proved by showing that the number of clans on the membrane drops 
to 1 at certain times. For small one would observe a cluster at these times. However the 
frequency of these events is proportional to A^~^, which indicates that polarity cannot occur 
in the large population limit — > oo, unless other mechanisms are present. 

In [19] we rigorously study this model under a different scaling of parameters. We multiply 
/cfb and /coff by A^, leaving kon and D unchanged. We keep track of the locations of the 
membrane molecules as well as their clan identities. A clan identity is a number in [0, 1] which 
is passed unaltered from the parent molecule to the offspring. A molecule that associates 
itself on the membrane is assigned a uniformly chosen clan identity in [0, 1]. At any time, the 
molecules on the membrane that have the same clan identity should have a common ancestor 
and hence they must belong to the same clan. Note that here the number of types (g) is 
equal to 1 and the population density is the same as the fraction of cell molecules that are 
on the membrane. Let E = E x [0, 1], where E is the membrane (sphere of radius R in M^). 
When the number of molecules is A^, the population dynamics is described by a AiN,a{E)- 
valued process : t > 0} as before. For any h G M+ let = A;fb(l — h), p{h) = kos 

and K,{h) = kon{^ — h). Let G V{E) be the uniform distribution on E and let the spatial 
migration operator B (see Section 2.1) be {D/2)A, where A denotes the Laplace-Beltrami 
operator on the sphere E. For any / : x [0, 1] — )■ M, A acts on / only as the function of 
the first coordinate. The operator {D/2)A is just the generator of the Brownian motion on 
E with diffusion rate D. With this notation one can verify that this model is a special case 
of the model in Section 2.3.4. The association mechanism is analogous to immigration while 
the feedback mechanism gives rise to the density regulation mechanism. Theorem 2.6 (see 
also Theorem 2.3 in [10]) shows that as A^ — j- oo we have /i^ =^ hcqi^ where 

^off 



kih 



and {iy{t) : t > 0} is a. P(i?)-valued Fleming- Viot process with generator 



p. rn , 

1=1 j^l l<l^k<m j^l,k 



+ 



for any F{y) = YYiLiifi^ ^) ^ ^o- It can be shown that this Fleming- Viot process is ergodic 
and has a unique stationary distribution in V(V{E)) (see Section 5 in [13] and Proposition 
2.5 in [19]). To study the evolution of the clan sizes we define a V{[0, l])-valued process 
Mt) : t > 0} by 

u,{t, S) = u{t, Ex S), 5 G i3([0, 1]). 

This is a Fleming- Viot process that describes the infinitely-many-neutral-alleles model (recall 
the discussion in Section 3.1). Therefore for any t > 0, we can write Ucit) = Si^i ca^xi, which 
means that Oj fraction of the population has clan identity Xi. At stationarity, the clan sizes 
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(arranged in descending order) are distributed according to the Poisson-Dirichlet distribution 
with parameter a = fcon/^fb- Properties of the Poisson-Dirichlet distribution (see [15]) tell 
us that for any small e > 0, there is a positive probability of the largest clan having size 
greater than (1 — e). Furthermore one can show that at stationarity the molecules in each 
clan are concentrated on a circular patch on the membrane. The square of the radius of this 
patch can be approximately computed as (see Theorem 2.7 in [19]) 



The last two assertions imply that if D is small in comparison to R^, then at stationarity 
there is a positive probability that most of the membrane molecules are in one clan and that 
clan is spread over a small area on the membrane. Due to ergodicity this event will occur 
infinitely often in any trajectory of the process {i^it) : t > 0}. Whenever this event happens 
we can expect the cell to be polarized. Therefore the limiting process exhibits recurring 
cell polarity. In [19] we discuss how the frequency of observing polarity depends on various 
model parameters. 

The above analysis shows that if the feedback mechanism is strong enough, it can counter 
spatial diffusion to generate cell polarity. However this conclusion is based on a highly sim- 
plified model. As mentioned in Section 1, most cells that exhibit polarity have complicated 
feedback circuits, with molecules of several types pulling each other on and off the membrane 
at various type-dependent rates. These different types of molecules may also have their own 
migration and association mechanisms. It would be interesting to know if the above analysis 
can be extended to general multi-type models for cell polarity. The results in this paper show 
that this can indeed be done as long as the feedback mechanism satisfies the assumptions in 
Section 2.2, and acts at a faster timescale than the association and migration mechanisms. 
In this case. Theorem 2.6 guarantees convergence to a Fleming- Viot process and this limiting 
process can then be analyzed in the same way as in [19]. This enables us to draw similar 
conclusions about the onset of cell polarity in this multi-type setting. 

4 Proofs 

4.1 Well-posedness of the martingale problems for Af^ 

Recall the definitions of the operators A^, Af, and Af from (2.22), (2.23), (2.24) and 
(2.26). In this section we prove Proposition 2.4 which says that the martingale problem 
for these operators is well-posed in the space D_x4<i^ (£;)[0, oo). Pick a / G {0, 1, 2, 3}. If we 

do not allow the dynamics under Aj^ to leave a compact set of Ai%^^{E), then Af^ can be 
viewed as a bounded perturbation of the migration operator (given by (2.12)). The well- 
posedness of the corresponding martingale problem is immediate from the well-posedness 
of the martingale problem for (see Chapter 4 in Ethier and Kurtz [12]). In our case, 
the dynamics under Af may exit any compact set of Ai%^{E). However we can still argue 
the well-posedness of the corresponding martingale problem by showing that this exit time 
tends to infinity as the compact set gets bigger and bigger in size. We now make these ideas 
precise. 



2D 
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Lemma 4.1 Fix a I e {0, 1,2,3} , N e N and n e V {M%^^{E)). For each k e N let 
{fj'kit) '■ t > 0} be a M.'^j^ J^E) -valued process with initial distribution tx. Define a stopping 
time 

r, = inf {t > : > k or > k} , (4.49) 

where H is the density map (2.2). Suppose that for each G N and F G Cq 

F{fik{t A n)) - F{fi,m - / Af'F{fik{s))ds 

Jo 

is a martingale. Then for any t > 

lim P {Tk <t) = 0. 

k—>oo 

Proof. Let {hk{t) = Hi^^kif)) : t > 0} be the density process corresponding to ^k and let 
Ck : ffi+ — > ffi+ be the function defined by 

c.W = (c„W.....c,,(ft)) = {l^ 'X'lwtee'* (^-^O) 

Pick any e G (0, 1). Since for each A; G N the distribution of //^(O) is tt there must exist 
a fee > such that 

P(||/ifc(0)||i > A;,) < e for all A; G N. 

If satisfies /c^ < ke then 

E((cfc(/ifc(0)),/ifc(0))) < K^{\\hk{Q)\\i < ke) + 2A;P(||/ifc(0)||i > A;,) < 3ek. (4.51) 

For each i G Q and /c G N, define a function F^'^ G Cq by F^{ij) = Ck,i{h)hi where 
h = H{fi). From (2.31) we know that for any fi G ^{E) 

AfFtif^) = B'^Ftif^) + iVR^^^/i) + Gf ^^/i). 

One can easily verify that B^_F;'=(/i) = 0, K^Fj'iij,) = Ck,i{h)ei{h), G^^'=(/i) = 0, G^^^(/i) = 
0, G^_F^/=(/i) = Cfc,i(/i)Ki(/i) and finally 

Gf^'=(/x) = CkAh) {^j^A^,h)i^,{dx) - l^dt{x,h)ft,{dx)^ . 

Note that for each i,j G Q, the functions are bounded, while the functions 6i and Kj 

satisfy (2.16) and (2.25). This implies that there exists a positive constant C (depending on 
and /) such that 

Af'F'^ifx) < C (1 + {Ckih), h)) for all ^ G M%^,iE). (4.52) 
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By the assumption stated in the statement of this lemma we can say that 



Jo 

is a martingale starting at 0. Taking expectations we get 

E [FtMt A r,))) = E [f^'MO))) + E Af" Ft{fids))ds^ • (4.53) 

Let F*-' : Mp{E) M be given by = '^i^zq Pf'dj,) = {ck{h),h). Then summing over 

i E Q in (4.53) and using (4.52) we arrive at 

From (4.51) and Gronwall's inequality, for k > k^/e we obtain 

E (F^(/ifc(t A Tfc))) < {3ke + Cqt) e^"*. 
Then by Markov's inequality 

EfF'=(/i,(tAr,,)) 



lim P (F''{fik{t A Tfc)) > A;) < lim 

A;->oo V / fc— >oo k 

^ (3A;e + Cgt) e^^* 
< lim 



Observe that 



fc— >oo k 
< 3ee^^*. (4.54) 



^in<t)=Fi\\hkitArk)\\, > k) 

= F(F\fikitATk)) > k) +F{\\hk{tATk)\\^ > 2k). 



For large k, the second probability on the right is because of the following reason. The 
process hk has jumps of size and hence the definition of (see (4.49)) implies that 
\\hk{t A Tk)\\^ <k + (l/N) < 2k. Therefore using (4.54) we get 

lim P(rfc < t) = lim P (P^fikit A r^)) > A;) < 3ee^''*. 

k—^oo fc— )-oo V / 

Letting e — t- proves the lemma. □ 

Proof. [Proof of Proposition 2.4] Fix a e N and a / G {0, 1, 2, 3}. The space M% ,^{E) is 
complete and separable and for each k E N the set 

U, = G M%^,{E) : ||iJ(/i)||i < k} (4.55) 
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is open with a compact closure in Jlip{E). Define an operator : V(Lk) = Cg — ;> 
B (Ml^aiE)) by 



L,F(/i) = B^F(/i) + (iVR^F(/i) + Gf^F(/i)) , 

for any F G Cq. The operator can be seen as a bounded perturbation of the operator B^. 
We argued in Section 2.1 that the martingale problem for is well-posed. From Theorem 
4.10.3 in Ethier and Kurtz [12], the martingale problem for is well-posed for each G N. 
Pick a vr G P (A^^^(i?)) and let {fik(t) : t > 0} be the unique solution to the martingale 
problem for (L^, tt). Define a stopping time by 

Tk = mi{t > : /ifc(t) ^ Uk or yUfc(t-) ^ f/J. 

Then for any F G 

- F(/ifc(0)) - / LkF{fik{s))ds 
Jo 

is a martingale. From (2.31) one can see that if fi E Uk then A[^i^(/x) = LfcF(/x). Using the 
optional sampling theorem we get that 

F(/iA:(t A Tfc)) - F(/ifc(0)) - / LfcF(/ifc(s))rfs 

^0 

= Fifikit A Tfe)) - F(/ifc(0)) - / A^F{fikis))ds 

Jo 

is a martingale. Lemma 4.1 ensures that for any t > 0, 

lim ¥{Tk <t) = 0. 

From Theorem 4.6.3 in Ethier and Kurtz [12] we can conclude that there exists a unique 
solution to the martingale problem for (A^,7r). □ 

4.2 Properties of the limiting process 

The goal of this section is to prove Proposition 2.5 which gives important insights into the 
limiting behaviour of the dynamics under the models described in Section 2.3. As mentioned 
in Section 2.4 our proof of Proposition 2.5 will rely on the work of Katzenberger [23] which 
studies how semimartingales behave when they are driven by a fast drift that forces them to 
stay on a stable manifold. Before we can use the framework in [23] we need to prove some 
preliminary results. We start by recalling a tightness condition for semimartingales. 

Condition 4.2 Let {Z^ : N E N} be a sequence of {J^^}- semimartingales with paths in 
-DKd[0, oo). Assume that for all T > 

sup \\Z^{t) - Z^{t-)\\ ^0 as N ^oo. (4.56) 

0<t<T 
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Moreover assume that each semimartingale can he written as 

Z''{t) = M^{t)+ f £'{s)ds 
Jo 

where is a square integrable {J^^^ } -martingale and is a {J^^} adapted process satis- 
fying ^ 

supe([M%+ f \\A^(s)\\ ds] <oc 



for each t >0, where [M ]t is the quadratic variation of the martingale M . 

Remark 4.3 // a sequence of semimartingales {Z^ : G N} satisfies Condition 4-2 then 
this sequence is tight in D^d[0, oo) in the Skorohod topology (see Corollary 2.3.3 in Joffe and 
Metivier [JI]) and any limit point Z is a semimartingale with continuous sample paths. 

Lemma 4.4 Pick a I E {0,1,2,3}, i E Q and f eVq. For each N eN, let : t > 0} 

be a solution of the martingale problem for Af. Define a real-valued process {Z^{t) : t > 0} 
by 



^''W=(/,-"fW)-(/,/if(0)) 



.N 



N I 



N 



J2 (3Ah''is)){f, /if (s)) - p.ih^is)){f, /if (s)) 



ds. 



where h'^lt) = H{n^{t)). Then Z^ is a semimartingale with respect to the filtration gener- 
ated by {l2^{t) : t > 0}. For any compact K C M!^ define 



\^{K) = inf <{ t > : /i^ (t-) iK or (t) 



N/ 



(4.57) 



where K denotes the interior of the set K. If Z^ is the semimartingale given by Z^[t) = 
Z^ {t A {K)) for t > 0, then the sequence of semimartingales {Z^ : N E N} satisfies 
Condition 4. 2. 

Proof. Let i E Q and / G be fixed. For each N eN and / G {0, 1, 2, 3} define a function 
af : W^iE) ^ M by 

«o^(/^) = {Bif,pi), 

O'lilA = {Bif,iXi) + ^ / bji{x,h)f{x)pj{dx)- / d'^{x,h)f{x)fii{dx), 
j&Q J^ J^ 



af (/i) = {B,f,^^,) + NY, f Pji{h)pf,{x) \ [ 

jGQ Lis 

and af (/i) = {Bif,pi) + Ki{h) / f{x)Qi{dx), 

Je 



{f{y)-f{x))K,{x,dy) 



/ij (dx) 
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where h = H{ii). Let C M%SE) be given by 

t/f = {/i G M^^SE) : < 2fc} . (4.58) 

Then for each /c G N and / G {0, 1, 2, 3} we have 

sup sup {jji) < oo. (4.59) 

iVGN ^gC/Af 

To see this note that for any iV G N, f/f C Uk := {yU G M^E) : ||i/(/i)||i < 2A;} and Uk is 
a compact subset of J^p{E). For / G {0, 1,3}, af^ is a continuous function which does not 
depend on TV and hence we get (4.59) simply by observing that 

sup sup (fi) < sup af^(/i) < oo. (4.60) 

Similarly if we define a continuous function a2 '■ Ai'p{E) — > M by 

a2(yu) = {Bif.jii) + ^l3jiih){Cjif,fij) 

then we also have 

sup sup a2(/i) < sup a2(/u) < oo. (4.61) 

Here Cjj's are the operators satisfying Assumption 2.2. This assumption also implies that 
sup sup |a^(yu) — a2(yu) I 

<sup sup V/3,.(/^) / (np%{x) [ ifiy) - fix))^f^{x,dy) - C.Jix)) fx,{dx) 

< sup sup ^f3ji{h)hj sup ( Npfi{x) / (/(y) - f{x))i^fi{x,dy) - Cjif{x) 

< OO. 

This bound along with (4.61) and the triangle inequality shows (4.59) for / = 2. 
Let Ck-.Rl^ M.I be given by (4.50). Define : M'^p{E) ^ M by 

FkifJ') = Ck,i{h){f,Hi), 
where h = H{fi). One can verify that for any fi G 

(B^ + Gf ) F,(/i) = (/i) (4.62) 
and R^'Fkifi) = ^ /3,,(/i) p,(/i)(/, /i,). (4.63) 
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Suppose : t > 0} solves the martingale problem for Af^ and {h^{t) = H{fi^{t)) : t > 

0} is the corresponding density process. Define another process {m^(t) : t > 0} by 



m^(t) = - / af'(/i^(s))rfs. (4.64) 
Jo 

For any /c G N let 

= inf{t > : \\h^{t-)\\i > k or \\h^{t)\\i > k}. (4.65) 

From Lemma 4.1 we can conclude that for any fixed N, the stopping times converge to 
oo a.s. as /c — )■ oo. Observe that belongs to the class Cq = P(Af^). Hence 

Ffc(/.^(t)) -F,(/x^(0)) - rAfF,(/x^(s))ds 

^0 

is a martingale. From (2.31) and the optional sampling theorem we get that 

<{t) = F,{f,^{t A rf )) - F,(/.^(0)) - / AfF,{f,^{s))ds (4.66) 
= Ffc(/i^(t A rf )) - Ffc(/i^(0)) - N / R^Ffc(/i^(s))rfs 
(B^ + GD F,.(/i^(s))rf. 



is also a martingale. If the set (0, t A t^] is non-empty then for any s G (0, t A r^], we have 
Ck{h^{s)) = Ig and therefore Ffc(/i^(s)) = (/, /if^(s)). If the set (0,t A r^] is empty then 
t At^ = and in this case (t) = 0. From (4.62) and (4.63) we see that for all t > 0, 

[t) = [t At^), where is defined by (4.64). But is a martingale and for a fixed 
N, ^ oo a.s. as /c — 7- oo. Therefore we can conclude that is local martingale and 
hence is a semimartingale. 

Let : M%{E) ^ R be given by F^{iJ,) = (-Ffc(/i))^. Note that for any /i G A^i.(E) 



AT -^7 AT Ar2 

Using this one can verify that if /i G f/^ and h = H{fi) then we have 



N (k'^PM - 2Ffc(/x)R^F,(/x)) = 5^/3,.(/^)(/^/x,) + p.(/^)(/', /^.), 
G^F'M-2FMG^hif^) 

= TtY. I b%ix,h)f{x)fi,{dx)+ I dt{x,h)f{x)fii{dx), 
G^F,2(/i)-2Ffc(/i)G^Ffe(y^) 
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and G^F-^ifx) - 2Ffc(^)Gf = l^,{h) j^f\x)ei{dx). 



Also for any /i G f/^ 

BV2(/i) - 2Ffc(/i)B^Ffc(/i) = - 2/5,/, /i,). 

Using (2.31) and the above expressions we can show in a manner similar to (4.59) that for 
any G N we have 

sup sup |Af F2(/x) - 2FMAf'FM\ < oo. (4.67) 

iVgN ^gC/Af 

The function F| is also in Cg = V[Af). Therefore if is the martingale given by (4.66) 
then 

Kit))' - {AfF^if^^is)) - 2h{^,^{s))ArF,{^,^{s))) ds 

is also a martingale. Therefore the expected quadratic variation of can be computed as 

E = E ((mf (t))2) = E (Ar^f (/i^(.)) - 2F,(/i^(.))ArF,(/i^(.))) rf. j . 

(4.68) 

For any fixed /c G N, the sequence of semimartingales {2'^(- A t^) : G N} satisfy (4.56) 
because the discontinuities of fi^ are of size proportional to 1/A^. From (4.64) we can see 
that the semimartingale Z^(- A t^) can be decomposed as 

Z%t A rf ) = m\t^ rf ) + / ' [^\s))ds = m^t) + / l^.^^.^af {^,'' {s))ds. 

Jo Jo 

For any < s < , /i^(s) G t/f . Using (4.68), (4.67) and (4.59) we can see that 

supEfK]i+ / l^,^,Nyaf'ifi''{s))ds) <oo. (4.69) 
Nm \ Jo J 

Therefore for any A; G N, the sequence of semimartingales {2'^(- A t^)} satisfies Condition 
4.2. For any compact set K C M^, there exists a k such that if C {/i G : < A;}. If 
the stopping time A;^ is defined by (4.57) then < a.s. where is given by (4.65). 
Hence it is immediate that if is the semimartingale defined by Z^{-) = Z^{- A A^(A')) 
then the sequence of semimartingales {Z^ : N G N} will also satisfy Condition 4.2. □ 

Proof. [Proof of Proposition 2.5] Let {J^l^} be the filtration generated by the process {/i^(t) : 
t > 0}. For any compact set K C Ueq let A^(A') be given by (4.57). In this proof a sequence 



33 



of {J-'/^}-semimartingales {Z^ : G N} with paths in D^d[0, oo) will be called well-behaved 
if for each compact K C Ueq, the sequence of semimartingales {Z'^{- A \^{K))} satisfies 
Condition 4.2. 

Using Lemma 4.4 with f = 1e, (where 1^; is as in (2.2)) for each i E Q we obtain a 
well-behaved M- valued semimartingale Zj^'^ such that 



(U,/xf(t)) = (U,/.f(0)) + iV 



J2 (3Ah''is)){lE, /if (.)) - P^ih^m^E, /if (S)) 



ds 



Recall the definition of the matrix A{h) from (2.13). The above expression is the same as 



J2PAh''is))hfis)-p,ih''{s))h^is) 



ds + Z^'\t) 



ds + Z^'\t) 



If we let Z^'^ to be the M'^-valued semimartingale given by 



Z^'i(t) = (Zf'^(t),...,Zf'i(t))fort>0, 
then Z^'^ is also a well-behaved semimartingale. The semimartingale h'^ satisfies 



(t) = (0) + N / A{h^{s))h^{s)ds + Z^'\t) 



h''{0) + N e{h^{s))ds + Z^'\t) 



(4.70) 



where the last equality holds due to definition (2.14). Now fix a / G Vq. From Lemma 4.4, 
for each i E Q there is a well-behaved semimartingale Z^'^ such that 



ieQ 



(/,/^f(t)) = (/,/^f(0)) + iV( j Y.f^dh''{s)){f,pf{s))ds- I p.{h''{s)){f,fif{s))ds 



(/,/if(0)) + iV 



5^A,,(/^^(.))(/,/if(s)) 



ds + Z^^f{t). 



Using the integration by parts formula for semimartingales, for each i,jEQ we can write 



C(t)(/,<(t)) = /^f(0)(/,/if(0))+iVl / h^s 



.N I 



N I 



.k&Q 



ds 



+ / (/,/^f(^)) 



J2A,{h''{s))h^{s) 
.keQ 



ds 1 + Z^ (t) 



(4.71) 
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where Zj^ is another well-behaved semimartingale given by 



ZS{t)= [ hf{s)Zf^\s)ds+ I {f,f,f{s))Z^'\s)ds + [Zf'\zf'^^ 
Jo Jo 

The last term in the above equation is the cross-variation term between Z^'^ and Z^'^ . Now 
for each i ^ Q define the semimartingale Y/^ by 



(4.72) 



The semimartingale Y/^ is just a linear combination of the semimartingales of the form 
hf{t){f,nf{t)). Using (4.71) we can write 



Y^t) = Y^iO) + N f Lf{s)ds + ZP(t), 
Jo 



N 2 

where Zj ' is a well-behaved semimartingale and for any t > 

Lfit)= E [A^k{h\t))h^{t){f,^,Ut)) + A,,{h''{t))hUt){f.^^nt)) 

j,k€Q 

-A,.(/^^(t))/.f (t)(/,/if (t)) - A.(/^^(t))/.f (t)(/,/if (t))] 



5^A,,(/^^(t))n^(t) + (/,/if(t)) 

-/zf(t) 1] A,.(/.^(t))(/,/.^(t)) 

.j,k€Q 



J2 A.^ih"^ {t))hUt) 

.j,k€Q 



Define a matrix G{h) for each h G M!!^ as follows 
G{h) = 



if /i = 0„ 



A(/i) 



^I^-^^UW) otherwise. 



Note that if h^{t) ^ Og then 

y,^(t) + /^f(t)(E,,Q(/,/if(t)) 

(/,^f(t)) = ^ and 



5^ A,,(/z^(t))(/,/if(t)) 
j,fceQ 



(4.73) 
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This allows us to write 



and hence from (4.73) 



N I 



ds + Z^'\t). 



(4.74) 



Let and Z^'"^ be the M"^ ^-valued semimartingales given by 



Y'\t) = {Y{^it), . . . , Y^^_,{t)) and Z^'\t) = (Zf '^(t), . . . , Z'Jltit)) for t > 0. 
For each h G let G(/i) G M(g — 1, g — 1) be the matrix defined by 

Gij{h) = Gij{h) - Giq{h) for all i, j G {1, . . . , g - 1}. 
Observe that EieQ^i^(^) = and hence Y^ [t) = -J2!ZlY-^{t). From (4.74) we get 



7N,2, 



Y^{t) = Y^{0) + N / G{h''{s))Y''{s)ds + Z''''{t). 



(4.75) 



From part (E) of Lemma A.l the matrix G{hcq) is stable, that is, all its eigenvalues have 
strictly negative real parts. 

We now define a x M'^^-'^-valued semimartingale by 

X^{t) = (/i^(t), y^(t)) = «(t), ...,h^{t),Yf{t),..., F^^i(t)) for t > 0. (4.76) 

From (4.70) and (4.75) we can see that satisfies 



X'' {t) = X'' (0) + N / F{X'' {s))ds + Z'' {t) 

Jo 

where Z^(t) = (Z^'^(t), Z^'^(t)) is a well-behaved semimartingale and F : 
M.'^i~^ is the function given by 

F{x) = {e{h),G{h)y) for x = {h,y) G x R'^-\ 



(4.77) 



(4.78) 



Let Xcq = (/ieq,Og_i). Thcu F(xeq) = 02g-i and the Jacobian matrix [JF(xcq)] G M(2g — 
l,2q — 1) has the block lower-triangular form 



C G(/leq) 



where C is some (g — 1) x g matrix in M(g — 1, g) and Og^g^i is the g x (g — 1) matrix of 
zeroes. We mentioned above that the matrix G{heq) is stable and part (B) of Assumption 
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2.1 says that the matrix [JO{hcq)] is also stable. Due to the block triangular form, the matrix 
[JF(xcq)] is stable as well. 

Pick any ho G f/cq and yo G M.'^^'^. Let h{t) = ipeiho^t) for all t > 0, where V'e is the flow 
defined in Section 2.2. The set Ueq C in ^/'g-invariant and hence h{t) G f/eq for all t > 0. 
Let be the unique solution of the initial value problem 

^ = G(Mi))?/, 1/(0) = yo. (4.79) 

Since the above differential equation is linear in the y variable, the solution y{t) is defined 
for all t > 0. Moreover h{t) — )■ /leq and G{h{t)) — )■ G{hcq) as t — j- oo. The matrix G{hcq) 
is stable and therefore y(t) — )■ Og_i as t — oo. Let U{xcq) = f/cq x R'^""'^. For each xq = 
{ho,yo) G f/(xeq) and t > 0, let ipF{xo,t) = {h(t),y(t)) with h(t) = ipe{ho,t) and ?/(t) being 
the solution of (4.79). 

The mapping ipp : U{xeq) x R+ — )• f/(xeq) is the flow of the vector field F on U{xcq). For 
all X G ?7(xeq) and t >0, i(:p satisfies 

i)F{x,t) = X + f F{tpF{x,s))ds. (4.80) 

From the discussion in the preceding paragraph we can conclude that 

lim ipp(x,t) = Xeq for all x G f/(Xeq). (4-81) 

We have assumed in this proposition that there is a compact set Kq C f/cq such that /i^(0) G 
-ffo a.s. for each iV G N. Note that any function f E T>o is bounded. The definition 
of the semimartingale guarantees that there is a compact set Ki C M'^~^ such that 
X^(0) = (/i^(0),F^(0)) G i^o X Ki C f/(xcq) a.s. for each N eN. 

Now consider the equation (4.77). For large values of A^, the semimartingale X'^ is driven 
by a large drift term of the form NF^X'^ {■)). The vector X IS 3j stable fixed point for this 
drift term and U{xeq) is its region of attraction. If we start in this region of attraction, 
then this drift is very forceful. It completely overwhelms the effect of the well-behaved 
semimartingale and drives X'^ to the stable fixed point Xcq. Moreover as gets large, 
the trajectories of X'^ start looking more and more like the trajectories of the deterministic 
flow ipp with time compressed by a factor of A^. These ideas are made precise in a much 
more general setting by Katzenberger [23]. We use Theorem 6.3 in [23] to deduce that for 
any T > 

sup \\X^{t)-iljF{X^{0),Nt)\\^0 8.sN ^oc, (4.82) 
te[o,T] 

where || ■ || is the standard Euclidean norm. From the definition of X^ and ■^/'f it is also clear 
that for any T > 

sup \\h^{t) -i)e{h^{0),Nt)\\^0 as N ^ oo. (4.83) 

tG[o,r] 

From now on, for any x G and e > 0, let B^{x) denote the open ball in centered 
at X with radius e. We have already argued that the Jacobian matrix of F at Xeq is stable. 
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From a simple linearization argument (see for example the proof of Theorem 3.7 in Khalil 
[24]) we can see that there exists a 5o > such that the open ball B'^^~^{xeq) C U{xcq) and 
for every 6 G (0, 6o) there exists a ?/^i?-invariant open set Ws whose closure Ws is contained 
in i?^^"^(xcq). From (4.81) we obtain 

Ko xKic\J{xe t/(xeq) : t) C Wg} . 

t>0 

As Ws is -i/'iT'-invariant, the open sets on the right are getting bigger and bigger as t increases. 
Compactness of Kq x Ki implies that there exists a > such that ipF{x,t) G Ws for all 
X G Kfj X Ki and t >ts. This immediately gives us 

sup sup Wi^riXyt) — XeqW < S 

and letting 5 — )■ we obtain 

limsup sup \\iPf{x, t) — XcqW = 0. (4.84) 

t— >oo x&KqxKi 

Now let tN he any sequence satisfying the conditions of this proposition. Then for any T > 



cq 



sup ||X^(t + ^Tv) — X, 
i6[0,T] 

< sup \\X^{t + tN) -lljF{X^{0),N{t + tN))\\+ sup \\l/jF{X^{0),Nt + NtN)-Xeq\ 

te[o,T] te[o,r] 

< sup \\X^{t + tN) -ipFiX^{0),N{t + tN))\\+ sup sup Wippix^t) - Xc 



te[0,T] t>NtNxeKoxKi 



where the second inequality is true because X^{0) e KqX Ki a.s. for each iV G M. From 
(4.82) and (4.84) we can see that as — t- oo 



sup \\X^{t + tN) -XeqW^O 

te[o,T] 



which of course implies that 



sup ||/i^(t + t7v) - /icq|| ^ and sup ||F^(t + tjv) || ^ 0. (4.85) 

t6[0,T] t€[0,T] 

This proves part (A) of the proposition since the norms || ■ || and || ■ ||i are equivalent in M''. 
From the definition of Y/' we can check that for each i,j G Q and t > 

The limits (4.85) immediately give us part (B) of the proposition for any / G Vq. As is 
dense in C{E), part (B) holds for any / G C{E). □ 

The following lemma will be useful in proving Theorem 2.6. 
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Lemma 4.5 Let the notation and assumptions be the same as in Proposition 2.5. Then 
there is a compact set K C f/eq such that for all T > 

lim P (h^{t) ^ K for any t e [0, Tl) = 0. 

Proof. Since the Jacobian matrix [J6{hcq)] is stable (part (B) of Assumption 2.1), by a 
linearization argument similar to the one referred in the proof above, we can find an e > 
such that the open ball B^{hcq) C f/gq and there exists a ?/^e-invariant open set such that 
its closure Ue C B^{hcq). One can argue as before that since Kq is a compact set, there 
exists a such that for all t >t^ and x E Kq, ipe{x, t) E 11^. If we define Kq as 

Ko = {he f/eq : Mh,Q C fJJ , (4.86) 

then it is a -i/'g-invariant compact set containing Kq. Because we have assumed that h^{0) G 
Kq a.s. for all G N, we must have that for all t > 0, tljg{h^ {0),t) G Kq a.s.. But Ucq is 
open in M.^, (see Section 2.2) and so there is a 7 > such that 



K = <x E Ucq : inf \\y — x\\ < 7 

is a compact subset of Ugq. Observe that if for some t > 0, h^(t) ^ K then we must have 
that \\h^{t) -tlJe{h^{0),Nt)\\ > 7. Therefore 

lim P (h^{t) i K for any t G [0,T]) < lim P | sup ||/i^(t) - i)e{h^ {^),Nt)\ > 7 ) . 

The limit on the right is due to (4.83) and this proves the lemma. □ 



4.3 Solution to a system of partial differential equations 

Recall the discussion at the end of Section 2.4. To prove Theorem 2.6 we require a function 
A that allows us to construct a V{E) valued process {u^ {t) : t > 0} (see (2.44)) whose 
dynamics is well-behaved as A^ approaches 00. The goal of this section is to guarantee that 
such a function A exists. 

Specifically, we need to show that for some open set Ucq C M'^ containing Ucq (given by 
(2.20)), we have a function A G C^(f/cq,IR^) which satisfies the following: 

A^{h)A{h) + [JA{h)] e{h) = Oq for all h G Ucq, (4.87) 
(A(/i), h) = 1 for all h G Ucq (4.88) 
and A{hcq) = Vcq. (4.89) 

Here Vcq is defined in (2.34) and [JA(h)] in equation (4.87) refers to the Jacobian matrix of 
A at h. The significance of the above relations will become clear in Section 4.4. 

The major difficulty in solving (4.87) arises in the neighbourhood of heq. This is because 
0{hcq) = Oq (part (A) of Assumption 2.1), which causes degeneracy in the system. However 
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the next proposition shows that by employing power series expansions we can get around 
this problem and find an analytic solution to (4.87) in a neighbourhood of /icq. We later 
construct an open set Ucq containing Ucq and extend the solution over the whole Ueq. We 
also show that this solution has all the properties we desire. 

Proposition 4.6 There exists an open set V containing /igq such that the equation (4.87) 
has an analytic solution A on V satisfying (4.89). 

Proof. We first transform the equation (4.87) into another equation that is easier to work 
with. Let Ai, . . . , Aq be the eigenvalues of the matrix [J9{heq)]- We will prove this proposition 
under the assumption that all these eigenvalues are real. We later remark how the proof 
changes when they take complex values. 

We know from part (B) of Assumption 2.1 that Aj < for each i & Q. Pick an eg G (0, 1) 
such that 

Xi < -4eo for all i G Q. (4.90) 

Let Ml G M(g, g) be the matrix representing the Jordan canonical form of [J6{hcq)]. Its 
diagonal is occupied by Ai, . . . , Ag, while its super-diagonal entries are either or L All the 
other entries are 0. Let Pi G M(g, q) be the invertible matrix such that Pi[J9{hcq)]Pi^ = Mi. 
Let P2 = Diag(l, eo, eg, ... , e^"^), P = P^^Pi and M = P^^M^Pi. Then 

P[je{h,^)]p-' = M, (4.91) 

and M is just the matrix Mi with each 1 on the super-diagonal replaced by eo. 

In this proof, will always denote the vector of zeroes in M.*^. Since h^q G (see part 
(A) of Lemma A.l) and the map x H- /icq + P~^x is continuous, we can find a ro > such 
that for any x G B'^^^O) we have /icq + P~^x G M^, where -8^0(0) is the open ball in M'^ with 
radius ro centered at 0. For all x G B'^^{0), let A{x) G M(g, q) and 6{x) G M'' be given by 

A{x) = A^(/icq + P'^x) and e{x) = P9{hq + p-^x). 
Suppose /3 : U — )■ is a function which is analytic in an open set U C -8^0(0) 

X eU 

A{x)l3{x) + [J/3{x)]e{x) = (4.92) 

along with 

m = v,q, (4.93) 

where v^q is given by (2.34). If C is the image of U under the map x h- >■ heq+P~^x, then 
V is an open set containing h^q and the function A : — )■ M"^ defined by A(/i) = /3(P(/i — /icq)) 
is an analytic solution to (4.87) satisfying (4.89). Hence to prove the proposition it suffices 
to show that equation (4.92) has a solution /3 in some neighbourhood of which satisfies 
(4.93). 
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We will be using the multi-index notation to write the power series in q variables. For 
any multi-index a = (ai, 02, . . . , dq) G Nq let |a| = ai + 02 + ■ ■ ■ + and a\ = ai\a2l ■ ■ ■ aq\- 
For two multi-indices u = {ui, . . . , z/g) G Nq and a = (ai, . . . , aq) G Ng, we say that u < a if 
Vi < Ui for alH = 1, . . . , g and we say that v < a ii v < a and v ^ a. \i v <a then 



V ) v\{a - v)\ ' 
For any vector x G M'' and multi- index a = (ai, 02, ... , aq) G Ng define 



_ ai 02 aq 



and the differential operator Da as 

Qai Qaq 



The operator acts component-wise on matrix and vector valued functions. 
Consider the function (3 given by the power series 



where 7q, G M3 is given by 



P{x) =v,^+Y, 7aa:", (4.94) 

|a|=l 



^ (4.95) 



This function (3 satisfies (4.93). To prove the proposition it suffices to show that the vectors 
7o can be suitably chosen such that /3 satisfies (4.92) and there exists a positive constant C 
such that 



Jloo < Cl°' for all a G N^. (4.96) 



The last condition ensures the absolute convergence of the power series (4.94) in a neigh- 
bourhood of 0. 

Since ^(0) = d{hcq) = 0, if we plug x = in (4.92) we obtain 

l(0)/3(0) = A^{hq)v,q = 0. 

This is satisfied because of the choice of feq (see (2.34)). 

Applying the operator Da to equation (4.92) and using the product rule for multi- 
derivatives we get 



= Da i^A{x)P{x)j + Da [[jP{x)]e{x] 

= E ( " ) (^("-)^(^)) + E ( " ) (DAJ(3{x)]) [D^a-.Ax] 
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i'<a ^ ' {u<a,\a-v\>l} 

On rearranging we obtain 



{!^<a,|a-!^| = l} ^ ^ 



!^<a ^ ^ {v<a,|a-!/|>l} 



(4.97) 



For any j G Q, let Cj G Nq be the multi-index (0, . . . , 0, 1, 0, . . . , 0), with the 1 at the j-th 
position. Observe that if |a — z/| = 1 then v = a — Cj for some j G Q. Therefore 



E o:j[J9{x)]kj (L'(„_e,.+efe)/3(a;)) . 



j,keQ 



Note that [J6'(0)] = M (see (4.91)). This matrix has the eigenvalues Ai,...,Ag on the 
diagonal and either or eo on the super-diagonal. For each j = 2, . . . , g let = eo if 
M(j_i)j = eo and ej = otherwise. Then for a; = we obtain 



E ( " ) {DAJm)]) (^(«-.)^(o)) 



{u<a,\a~u\=l} 

Q 

= ^aj^jDaPiO) + E"iei^{"-ei+e,-i)/3(0). 
jeQ j=2 

Note that ^(0) = P6{hcq) = and for each a G Nq, 7„ is given by (4.95). We plug x = in 
(4.97) and divide by a\ to get 

q 

A{0)-fa + E "i^i^a + E ^"i-i + l)ej7(a-e,+e,_i) = (4.98) 

JGQ j=2,aj>0 

(a — z/j! z/! (a — uy. 

i'<a ^ ' {v<OL.\a~v\>\) ^ ' 
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The second term can be simplified as 

^ z/! (a - u)\ ^ ^ [oL-v)\ v\ 

{v<a\a~v\>\) ^ ' {v<(x\a-v\>\\ i^iQ ^ ' 



(a -I/)! 



Therefore we can write Yn as 



^" = -ES E E'^^ + l^l ^7(.+e,). 4.99 

u<ct ^ ' {v<OL.\a~v\>\) ^ ' 

For each /c G N, let Sk be the set of multi-indices given by 5^ = {a G Nq : |a| = A;}. The 
number of elements in is 

' A; + g - 1 
g-1 

We order the multi-indices in as follows. We say that v ^ a \i and only if XIigq — 
J2ieQ ^^i- Q^'^(l)! • • • ) (^^{^k) be all the elements of 5^ listed in the order given by ^. 

Let the matrix S*^^'^ G M(gsfc, qsk) be a block matrix composed of s\ blocks of size q x q. 
For each i,j G {1,2,..., s^} the block starting at row q{i — 1) -|- 1 and column q{j — 1) + 1 
of matrix H*^^^ is occupied by the matrix Lij G M(g, q) defined as follows. If i = j then 

La = A{0) + (^i^Q af{i)\i^ Iq. If i and j are such that a'^(j) = a''{i) — ei + ei-i for some 

/ G {2, . . . , g} then Lij = e;(af„^(z) + 1)/^. For every other i and j, Lj^ is just a matrix of 
zeroes. The matrix 3^^^^ is lower block-triangular and its determinant is given by 



Sk 



Det (SW) = n^et m + 5^af(^)A, ) /. 
i=i \ \leQ 



q 



The eigenvalues Ai,...,Ag satisfy (4.90). Since all the eigenvalues of the matrix A{0) = 
A^lheq) have non-positive real parts (see part (B) of Lemma A.l), the above determinant is 
non-zero. Hence the matrix H'^'^^ is invertible. 

Let X^''^ and Y^''^ be the vectors in M.''^'' given by 

X''''^ = (7afc(l); 7a'=(2); • • • ; lak(sk)) ^ud F^''^ = (^^^(i), yQ,fc(2), • • • , ^afc(sfe)) • 

Using (4.98) we obtain the following linear system 
and since the matrix S*^^) is invertible 

xik) = [s^j-iyw. (4.100) 
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Note that Y^''^ only depends on {70, : a E Si for / G {0,1, . . . , k — 1}}. Hence for each k eN 
we can solve for the whole set {7^ : a E S^} using (4.100). Doing this iteratively for each k 
we can solve for 7q, for all a E Nq. The function /3 given by (4.94) with this choice of ■jaS 
will solve (4.92) in a neighbourhood of if we can show that (4.96) holds for some C > 0. 
Showing this will be our next task. 

Any entry on the diagonal of E^''^ has the form Ajj(O) 4- XIjgQ ^j'^j ^'^^ some a E Sk and 
i E Q. Observe that A(0) = A^(heq) and this matrix only has non-positive entries on its 
diagonal (see (2.13)). From (4.90), for a G 5*^ we obtain the estimate 



jeQ 



> 4eofc. 



(4.101) 



For each row of HW, the sum of the absolute values of the non-diagonal entries is bounded 
above by 



max 



5^ 4(0) + J]q(«z-i + 1) 



1=2 



< max 



+ eo{k + q). (4.102) 



Hence from (4.101) and (4.102) we can conclude that there exists a. Kq E^ such that for all 
k > Kq the matrix S'^*^-' is strictly diagonally dominant and we have 



mm 

l<l<qsk 



r=l,rj^l 



lr\ 



> keo- 



Theorem 1 in Varah [36] shows that for all k > Kq 



1 



(4.103) 



Part (D) of Assumption 2.1 says that for each i,j E Q, the functions pi and (3ij are 
analytic in a neighbourhood of heq. This implies that there is a neighbourhood ?7 of such 
that the M(g, g)-valued function A and the M'^-valued function 6 are analytic component- wise 
on U. Therefore there is a constant Cq such that 



\\D^A{0)\\^ < and \\dJ{0)\\^ < C^^'a! for all a E Ng. 

We can assume that Cq > 1. Choose a 5 > satisfying 



6 < 



Coq{q + l)2'?+3 



(4.104) 



(4.105) 



and define C = Co/6. We will prove (4.96) by induction. Let k > Kq and suppose that C is 
large enough to satisfy 



(4.106) 



44 



for all / G {1, 2, . . . , — 1} and v G 5"^. To prove (4.96) we need to show that llTall^o < 
for all a G 5*^. This is equivalent to showing that ||X*^''''' ||^ < . From (4.100) and (4.103) 
we have 

II lloo — II L' ' J llooll lloo — /u6o 1 1 oo ' 

Hence to prove (4.96) it suffices to show that 

||r(^-)|L = max||F„|L<A;eoC"=. (4.107) 

From (4.99), (4.104) and (4.106), for any a G 5"^ we get 
But C = Cq/5 and |a| = /c. Hence 

W<a {u<a,\a—ii\>l} 

< j + 2kC ^ ^l"-'^! I . (4.108) 

\v<OL {i^<a,|o-!^|>l} / 



Note that 



1 + ^ = ^ = 5^ 5^ ■ ■ ■ 5^ n ^^"^"'^^ 

j/<a v<a ui=0 U2=0 Uq=0 i=l 



1 



S 



n i^t"" 



i=l,ai>0 



6 



and this shows that 

^l^-H < (1 _ _ 1^ (4.109) 



where n{a) be the number of non-zero coordinates of a. Similarly 

{u<a,\a-u\>l} v<ot 
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Since 5 G (0, 1/2) and n{a) < q, by Taylor's theorem we see that 

(1 - 1 < g2^+i5 and ^ ^l""'^! - 1 - n(a)5 < g(g + l)2'?+^52_ 

Using these estimates, (4.109), (4.110) and (4.108) we get 

\\Ya\\^ < (g2'^+i5 + 2kCq{q + l)2''+^6^) = CHk {C^q{q + 1)2^+^) _ 

But 5 satisfies (4.105) which shows (4.107) and completes the proof of the proposition. 

At the beginning of the proof, we had assumed that the eigenvalues Ai, . . . , of the 
matrix [J9{hcq)\ are all real-valued. If that is not true then the invertible matrix P that 
appears in (4.91) has complex entries. Let C be the field of complex numbers. Define a map 
: M"^ — 7- by 0(/i) = P{h — /icq)- The image of this map, denoted by (j){W^), sits as a 
g-dimensional real vector space in C^. The map (f) is an infinitely differentiable isomorphism 
between and </)(]R''') and using this we can define derivatives of real-valued functions over 
0(M'^). As above, we can obtain an analytic solution /3 of (4.92) satisfying (4.93), defined 
on some open set U in (f){M3) containing 0. On V = 0~^(f/), the function A defined by 
A(/i) = /9(0(/i)) will then be an analytic solution to (4.87) satisfying (4.89). □ 

The above proposition provides us with an analytic solution to (4.87) in a neighbourhood 
of heq. Our next task is to extend it to a solution in C^([/eq, 1^*) where Ueq is an open set in 
containing Ucq- 

Recall from Section 2.2 that for all i,j G Q, Pij,Pi are functions in C^(]R!^, IR+). Let 
O C M"? be the open set containing M!?^ defined by 

O = {h eM.'^ : hi > -1 ioY aWi = 1, . . . ,q} . 

Then we can extend the functions Pij, pi to functions Aj, Pi G C^(]R^, ]R-|_) such that (3ij{h) = 
and f)i{h) = for all h ^ O. Moreover since each Pij is bounded, we can make sure that its 
extension is also bounded. For each /i G let A{h) G M(g, q) be the matrix defined by 
(2.13) with Pij,pi replaced by Pij,pi. Also let 6 e C'^{W,W) be the function given by 

e{h) = A{h)h for h G R'^. (4.111) 

Corresponding to 6 we can define the flow map tp G C'^{W x M+, O) as the unique solution 
to the equation analogous to (2.18), with 9 replaced by 9. Define the region of attraction of 
the fixed point /icq as 

Ueq = \heO: lim i^{h,t) = kA . 
Then Ucq is an open set in (see Lemma 3.2 in [2 !]) containing Ucq- 

Proposition 4.7 There exists a solution A G C'^{Ucq,M.%) of (4.87) satisfying (4.88) and 
(4.89). 
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Proof. Suppose that U C Ucq is any T/^-invariant open set and the function A G C^(f/, M'') 
satisfies (4.87) and (4.89). We first show that this function automatically satisfies (4.88) on 
U. Using (4.87) and the -j/^-invariance of U we get 

j^HHh, t)) = [JAiiih, meii^ih, t)) = -A^iiih, t))A(^(/i, t)). (4.112) 
Observe that 

j^mh, t), A{4^{h, m = j^^ih, t),A{4^{h, t))^ + (i^ih, t), J t)) 

9(ip{h,t)^ ,A(V^(/i,t))) - [i^{h,t),A^{ip{h,t))A{^{h,t)) 
9(ip{h,t)),A{ip{h,t))) - (A{ip{h,t))ip{h,t),A{ip{h,t)) 



= 0, 

where the last equality holds due to (4.111). This shows that for any fixed h & U the 
function {ijj{h,t), A{ilj{h,t))) is a constant function of time. Therefore (4.89) implies that for 
any h & U 

{h,A{h)) = \im{^{h,t),A{^{h,t))) = (/leq,A(/leq)) = 1. 
f— ^oo 

This proves that A satisfies (4.88) on U. For any h & U and < t < let ^{h,t,tQ) be the 
matrix defined in Lemma A. 2. Since A satisfies (4.112) we must have 

A{^{h, t)) = $(/i, t, to)A(V^(/i, to)). (4.113) 

From Proposition 4.6 we know that on some open set C M'' containing h^q we can find 
a solution A G C^{V,R'i) that satisfies (4.87) along with (4.89). Since v^q G (that is, it 
is positive component-wise) and A is a continuous function, by shrinking V if necessary, we 
can ensure that the image of V under A lies in M^. Since V is open, there exists a r G (0, 1) 
such that B^{heq) C V, where B^{heq) is the open ball in centered at h^q with radius 
r. As in the proof of Lemma 4.5, we can find a -^-invariant open set W C B'^{hcq) which 
contains heq- 

For each n G N define an open set 

Or, = {heU,q:^{h,n)cW}. 

Each On is ■i/'-invariant. Furthermore W C Oi C O2 ■ ■ ■ and IJ^^i = Ucq- Define Xn{h,t) 
for each h G On and t G [0, n) by 

Xn{h, t) = $(/i, t, n)A{^{h, n)). (4.114) 

Observe that i/j{h,n) G W C V and so A{ip{h,n)) is well-defined and also A{ip{h,n)) G 
Wi. Part (C) of Lemma A.2 shows that A„(/i,t) G M«. Since A G C^iy.Wi), $(-,-,n) G 



47 



C2(f/,q X [0,n],MiR(g,g)) (see Lemma A.2) and e C^{W x R+,W') we must have that 
A„ G C^iOn X [0,n),M9). Note that A satisfies (4.113) for all h e W and so for < t < to 

A(^(/i, t)) = $(/i, t, to)A(V^(/^, to)). (4.115) 
Therefore ii h & W, then for any n G N and t > we have 

A„(/i,t) = A(z^(/i,t)). (4.116) 

Using parts (A) and (B) of Lemma A.2, (4.115) and the semigroup property of (similar to 
(2.19)) we can also see that for any h E On 



Xniil^iK t), 0) = $(V^(/i, t), 0, n)A [ilj{ij{h, t),n) 
= $(/i, t,n + t)A {^{h, n + t) j 

= $(/i, t, n)$(/i, n,n + t)K{4){h, n + t)) 
= ^{h,t,n)A{ij{h,n)) 
= K{h,t). 

Let h G 0„ and m > n. Then t/j{h,n) G IV. From part (A) of Lemma A.2 and (4.115) we 
can deduce that for any t G [0,n) 

Xm{h,t) = ^{h,t,m)A{ilj{h,m)) = ^{h,t,n)^{h,n,m)A{ilj{h,m)) 

= ^{h,t,n)A{i){h,n)) 
= K{h,t). 

Hence if we define the map A : Ucq x ]R_|_ — )• by 

X{h,t) = \n{h,t) if (/i,t) G 0„ X [0,n), 
then A is a well-defined function in C^(f/cq x M+jM^) which satisfies 

A(/i, t) = A(^^(/i, t), 0) for all (/i, t) G t/cq x M+. (4.117) 
From (4.114) and the definition of the matrix $ we can see that 

'^^^ = -A'm^))\{h^). (4.118) 

Define A : f/eq ^ M« by 

A{h) = X{h,0). 

Then this map is in C'^{Ucq,'^'i) and (4.117) implies that for any (/i, t) G t/. 



eq 



X 
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Using (4.118) we obtain 

[jK{^{h,t))]e{^{h,t)) = -^{i,{h,t))\{h,t) = -^{i,{h,t))K{^{h,t)). 

If we set t = then we see that A is a solution to (4.87). Since A satisfies (4.89), equation 
(4.116) imphes that A will also satisfy it. We have already shown that such a solution of 
(4.87) will automatically satisfy (4.88) for all h G t/eq- This completes the proof of the 
proposition. □ 



4.4 Fleming- Viot convergence 

In this section we will finally prove the main result of our paper, which is Theorem 2.6. Let 
{/x^(t) : t > 0} be a A^^^(i?)-valued process with generator A^^ for some / G {0, 1, 2, 3}. As 
outlined at the end of Section 2.4, we first extract a ■p(-E)-valued process {v^ {t) : t > 0} from 
the process {fi^{t) : t > 0}. This step requires a solution A of (4.87) whose existence was 
shown in Section 4.3. We then show that as — oo, we have ^ v where {^{t) : t > 0} is 
an appropriately defined Fleming- Viot process. This convergence and Proposition 2.5 prove 
Theorem 2.6. Before we proceed we need some preliminary results. 
For any set A C define 

M%{E : A) = {/i G Ml{E) : iJ(/i) G A] . 

Note that if A is a compact set then the set Ai'^p{E : A) is also compact. 

Recall the definition of the set f/eq from (2.20). Let {F/v : A^ G N} be a sequence of real- 
valued functions on A4p{E : Ucq). We will say that this sequence belongs to class o{N~'^) 
for some m G Nq, if and only if for each compact K C Ucq we have 

limsup sup N"" \FN{fi)\ = 0. 

For two such sequences {F^ : A^ G N} and {Gn : A^ G N}, we say that -F/v(yu) = GAr(/i) + 
o(A^-™) for all fi G M%{E : t/cq) if and only if the sequence of functions {{Fn-Gn) : A^ G N} 
is in the class o{N~''^). 

From now on let A G C2(f/cq, M^) be a function that satisfies (4.87), (4.88) and (4.89) on 
some open set t/eq C containing [/eq. Such a function exists by Proposition 4.7. Define a 
continuous map F : M''p{E : U^q) V{E) by 

F(^) = (4.119) 

where the measure v is given by 

u[S) = J2Hh)fiiiS) for any S G B{E), (4.120) 

ieQ 

with h = H{fi) being the density vector corresponding to fi. Note that for each h G Ueq, 
A{h) is a vector which is positive in each component and hence i^iS) > for all S G B{E). 
Since the function A satisfies (4.88) we have 

u{E) = J2Hh)f^i{E) = Y,Hh)h^ = 1. 
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This shows that u is a. probabihty measure on E. 

Let T' be the class of functions in C {J^p{E : Ucq)) given by 

r = = {h,{f, fi,) - h,{f, fi,)) L{fi) : ih, ...,h,) = H{fi), (4.121) 

/ G CiE) ,LeC (MUE : f/cq)) and t,j E Q} . 

Let T be the smallest algebra of functions in C {M'^p{E : Ucq)) containing T'. Observe that 
if G G C {Aip{E : Ueq)) and L G T, then the product GL is in T. Given two functions 
Gi, G2 G C (M'piE : f/eq)) we say that G2(/u) = G'i(/i) + T for all fi G M^E : f/cq) if and 
only if the function {G2 — Gi) is in the class T. 

Let E G G(V{E)) be a function in the class Cq defined by (2.4). Then E has the form 

m 

F{u) = l[{f,,u), (4.122) 

where /i, . . . , /m G Pq- Corresponding to E, define the functions Ei, Eik G Co for all distinct 
l,k eQhj 

m m 

Fi{u)= n (/i'^) ^^fcW= n (/^•''^)- (4.123) 

Using any F G Co we construct a function F G Cq as follows. We first extend the definition 
of A to the whole of W by letting K{h) = Og for all h ^ Ucq. If E has the form (4.122) then 
consider the function E : Ai'^p^E) M given by 

E{f^) = f[[J2A.ih){f,,fiA, (4.124) 

j=l \i€Q J 

where h = H{fi). Due to (4.88), the function E is in the class Cq defined by (2.11). The next 
result demonstrates how the action of various operators on functions of the form (4. 124) can 
be approximated. 

Proposition 4.8 Let E E Cq have the form (4.122). Corresponding to E let E E Cq have 
the form (4.124) and for distinct l,k E Q let Ei,Eik be given by (4.123). Then for all 
/i G A^^(£' : Ucq) with h = H{fi) and v = r(/i) we have the following. 

(A) Let be the operator given by (2.21). Then 

iVR^F(/i)= Yl lih)i{fifk,i^)-{fi,u){fk,u)))Eikiu) + T + oil), (4.125) 

l<l^k<m 



where 

^ " - " - (4.126) 



7(ft) = I 



J2 P,Ah)h,{hj(K)y- + pi(K)h,{Mh)Y 
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(B) Let B be the operator given by (2.12). Then 

BV(/i) = f2 ( J]A,(/i)/i.(S,/;,z/) I Fz(z/) + T + o(l). (4.127) 

1=1 \ieQ J 

(C) Let Gi be the operator given by (2.28). Then 

m 

Gf%) = Y.[{{b^{;h)fi{-),u) - {b^{.,h),u){fi,u)) (4.128) 
1=1 

+ {{d%; h), u) {fu v) - {d^i; h)fi{-), 1^))] Fi{u) + T + 0(1), 
where for any x E E and h G Mi 



h\x, h)=Y^ btjix, h)A.j{h)hi and d'{x, h) = J2 ^tix, h)Ki{h)hi. (4.129) 

i,j£Q i£Q 

(D) Let G2 be the operator given by (2.29). Then 

G^^(/i) = E ( E ) H^) + ^ + o{l), (4.130) 

1=1 \i,jeQ J 

where the operators Cij are as in Assumption 2.2. 

(E) Let be the operator given by (2.30). Then 

GfF(/i) = f](5^«:,(/i)A,(/i) [ ifiix) - {fi,u))Q,{dx)] Eiiu) + T + 0(1). 

1=1 \i&Q J 

(4.131) 

Proof. For any j G Q, let Cj be the vector in of the form Cj = (0, . . . , 0, 1, 0, . . . , 0) with 
the 1 at the j-th position. Since A G C^([/eq,ffi+) and f/eq is an open set containing [/eq, if 
h G t/eq, then using Taylor's theorem we can write 

Ki\h± — e,- = KAh) ± + '\ ' + oiN^^) 

"\ N ^ ) ^ N dhj 2iV2 Qi^2 ^ V ; 

for any z, j G Q. But then for any /i G Aip{E : f/gq) and x E E 



^A,(A)(/,.ft>±l(^ 



m 



n 

(=1 
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1 /lv:^9AK'!),, , aA,(ft),, A 

\ ieQ ^ ■> I 

^ m ^ m 

«=1 (=1 
+ ^ xi{l^,x)xiil^,x)Fik{u)+o{N-^), 



l<l^k<m 

where Fi,Fik are as in (4.123) and X/ (z^, a^), 0; (/^' ^) given by 

and 0K/^,^) = ^E^||^(/''/^^) + ^^/'(^)- (4.132) 

On rearranging we obtain 

-/1\- l"" l"" 

FU±-6i]- F{f,) =±j.J2^i{fi^ ^)^'(^) + ]^ E ^)^'(^) (4.133) 
^ 1=1 1=1 

+ ^ E x/(/i,a;)xi(/^,a:)F,,(i.)+o(iV-2). 

I<i5^fc<m, 

Therefore for any /i G M.p[E : Ucq) 



m / 

= E PM{xi{^^,■).^i^) - EpK^)(xi(/i,-),/i.) ) 

1=1 \i,jeQ ieQ 



+ E E P^Ah){^i{^,,.),^^,) + J2p^ih)m^i,■),^i.) ] F^i 

1=1 \i,jeQ i£Q 
l<l^k<m \i,jeQ 

+ E'^^(^)(^'('"' ■)^'('"' ■)""^) ) ^''^•(^) + ^(1)- 



But note that 



E (^iji^)(xi{f^,-),f^i) + 'Ypi{h){xl{f^,-),f^i) 
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,(/i)A,(/^) 



where the matrix A{h) and the vector 6{h) are defined by (2.13) and (2.14). Since the 
function A satisfies (4.87), the expression on the right is just 0. Hence the formula for 
iVR^F(/i) simphfies to 



= E ( E A.(/^)(0K/i,-),/^.) + ) ^'('^) (4.134) 

1=1 \i,jeQ ieQ / 

l<l^k<m \i,j£Q ieQ / 



Equation (4.88) says that for all h E 



eq 



J2hMh) = 1. 



Pick a j E Q. Differentiating the above equation with respect to hj we get 
and differentiating again with respect to hj we obtain 



Recall that for any ji G Aip^E : Ucq), v = r(/i) is given by (4.120). Using (4.88) one can 
verify that for any / G C{E) and i E Q 

if, ^>^ = h,{f,u) + Y,{h,{f,ii,)- h (/, h)) Hh) ■ (4.137) 

But the second term on the right is a function in the class T. Hence for all /i G AipiE : Ucq) 

{f,^i^)=h,{f,u) + T. (4.138) 

From the definitions of xj and (j)^ (see (4.132)) it is immediate that for any j E Q, I E 
{1, . . . , m} and x E E we have the following relations. For all yU. G A^_f(£' : Ucq) 



\ieQ ^ 



ibi 

Q 
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(//, ^) + T 



and ^Kfi,x) = ^^Mx)+(\j2h 



dAKh) 



Q J 



Using (4.135) and (4.136) we obtain 



and (pj (/i, x) 



dKj{h) 
dhi 



{fi{x)-{fuv)) + T. 



(4.139) 
(4.140) 



Recall that the class T is invariant under multiplication by functions in C{M.f{E : f/cq))- It 
can be checked that for any i,jEQ and /, /c G {1, . . . , m} 

ixiif^, Oxil/i, ■), /i.) = h^{A,{h)y {{fj,, v) - {fu v){h, v)) + T 

and the function /i H- (0^ (//,■), yU.j) belongs to class T. Substituting these two relations in 
(4.134) proves part (A) of this proposition. 

Recall the definition of the operator from Section 2.1. If G{v) = Y[]=i{9j^ ^ 
then one can verify (see Section 2.2 in [7]) that there is a constant c (depending on / and 
gi, ...,gi) such that 



sup sup I n 



I I 
j=l k=l,k^j 



< C. 



From the definition of the operator and (4.138) it is immediate that 

B^%) = Ul^^mB^fu^^^)\ Fi{v) + o{l) 

1=1 \ieQ / 

= E ( E ^■'(^)^^^^^^'' ^) ) ^'('^) + ^ + ^(1)- 

1=1 \ieQ J 

This proves part (B) of the proposition. 
Using (4.133), (4.139) and (4.138) we get 

Gf F(^) = N(j2f ^U^^ ^) (/^ + ^'^0 " ^^^^) ^'^'^''^ 
+ dtix, h) (f L - - F{^)\ ^,{d, 

m / 

- 5^ (rf|(-, h)xl{fi, ■),f^^)] Fi{u) + T + o(l) 
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1=1 \i,jeQ 

- J2 Mh) {dti; h) (/,(■) - u)) , /i,) ) FKi^) + T + o(l) 

m / 

= E E ^■(^)^^ (^^.(■' ^) (/'(•) - ^)) ' ^> 

- J2 Mh)h. {dti; h) (/,(■) - {fu u)) Fi{v) + T + o(l) 

m. 

= E(((^^(-'^)/'(-)'^)-(^^(-'^)''^)(/'''^)) 

Z=l 

+ /i), z/)(/,, z/) - z^))) F^z^) + T + 0(1), 

where the functions 6^ and are defined in (4.129). This proves part (C). 
Observe that Assumption 2.2 imphes that for any f E Vq and i,jEQ 



N J^p^^ix) /J/(y)-/(x))^,^}(x,rfy)/i,(rfx) = (a,/,/i.)+o(l) 



and 



Piji^) / {f{y)-fix))^^^{x,dy)fi,{dx) = o{l). 
Using (4.133) , (4.132) and (4.138) we obtain 

= N' (5,,{h) l^pf^ix) (^F (^/i + j^6i^ - F (^/i + j^6i}j^ ^i}(x, dy)^ f,,idx) 
= E E N{h) j^4{x) ^ (x^(/i,y) - xK/^,^)) ^5(x,rf|/)/i.(rfx) j FKz/) + 0(1) 



=1 \ jjeQ 



^(iV^A,(/^)A,(/^) / vU^)\l 

1=1 \ i,ieQ 



(/z(2/)-/z(a;))C(a;,rfy) 



/Xi(dx) FKz/)+o(l) 



= E E N{h)kAhm,fu^Ji^) Fi{u) + 0(1) 
1=1 \ijeQ / 

= E ( E A,(/^)A,(/^)/i.(C.,/z, u)] Fi{v) + T + 0(1). 

Z=l VijgQ / 

This proves part (D). 
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Again using (4.133), (4.139) and (4.138) we see that 



1=1 \ieQ ''^ J 

f;(^«:,(/.)A,(/.) / 




This proves part (E) and completes the proof of this proposition. □ 

Remark 4.9 For any {h,ij) G f/eq x V{E) let C^{h,u),CFih,'^),CF'\h,u)XF'^ih,u) and 
Cp'^ih, v) he the first terms that appear on the right hand side of equations (4.125), (4.127), 
(4.128), (4.130) and (4.131) respectively. Let C?'°{h,Jy) = for all {h,u) G t/eq x V{E). 
Recall the definitions of the operators A; fori G {0,1,2,3} from (2.39), (2.40), (2.41) and 
(2.42). One can verify that for each I E {0, 1, 2, 3}, F e Cq and v e V{E) 

AiF{v) = C|(/leq, U) + C|(/leq, ^) + Cf\K^. ^) ■ (4-141) 

To check this relation observe that A satisfies (4.89). Furthermore for each i G Q, Oiihcq) = 
"^ji^Q (^ji{heq)hcqj — Pi (/^-cq) ^cq,i = 0, wMch shows that thc valuc of the function 7 (given by 
(4.126)j at /icq is equal to the constant '^smpi (given by (2.36)). 

We now prove the main theorem of the paper. 

Proof.proof of Theorem 2.6] Fix a / G {0,1,2,3} and let {/i^(t) : t > 0} be a solution 
to the martingale problem for A^. Let {/i^(t) = H (/i^(t)) : t > 0} be the corresponding 
density process. Since /i^(0) =^ /i(0) as — )■ 00 and if(/i(0)) G Ucq a.s. we must have that 
h^{0) =^ h{0) and h{0) G Ucq a.s. It suffices to prove the theorem under the assumption 
that for all A^ G N, /i^(0) G Kq a.s. for some compact Kq C f/eq- 

By Lemma 4.5, we can find a bigger compact set K C Ucq containing Kq such that if we 
define the stopping time by 

= inf {t > : h^{t) ^ K} 

then 0"^ ^ 00 as A^ — )■ 00. 

Let F G Co be a function of the form (4.122) and let F G Cq have the form (4.124). From 
(2.31) and Proposition 4.8 we can conclude that for all /i G M.p[E : Ucq) 

Af^F(/i) = C|(/i, 1^) + Cf(/^, ^) + CF'\h, + + 0(1), (4.142) 

where h = H{fi), v G r(/i) and the continuous functions Cf, Cf, Cf'°' Cf'^; Cf'^; Cf'^ from 
f^cq X 'P(-E') to M are defined in Remark 4.9. 
The relation (4.142) implies that 

sup sup |AfF(/i)| < 00. (4.143) 
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The function F belongs to the domain of Af and the process {fi^(t) : t > 0} is a solution 
to the martingale problem for A^^. Hence 



F(/i^(t))-F^^(/i^^(0)) 



is a martingale and using the optional sampling theorem we see that 



F{fi'' {s))ds 



is also a martingale. Note that for all t > 0, fi^{t A o"^) is in the set M.f{E : f/cq). Define a 
P(£')-valued process by 



\t) = r A a^)) for t > 0. 



(4.144) 



The form of the functions F and F shows that F(fi^(t A o"^)) = F(iy^(t)) for any t > 0. 
Hence the martingale rrip can be rewritten as 



m^(t) = F(z/^(t))-F(z/^(0)) 



N/ 



tAa" 



Af'F(/i^(s))rfs. 



(4.145) 



The linear span of functions in the class Co is a dense sub-algebra of C {V{E)) and for every 
-F G Co we have the martingale relation (4.145). Theorems 3.9.1 and 3.9.4 in Ethier and Kurtz 
[12] along with the estimate (4.143) imply that the sequence of processes {u^ : G N} is 
tight in the space Dp(£;)[0, oo). 

Let tN be a sequence satisfying the conditions of Theorem 2.6. Pick a T > 0. It is easy 
to see that 



sup 

i6[0,T] 



AfF{^^''{s))ds - 



AfF{i,''{s))ds 



as cx). 



From parts (A) and (B) of Proposition 2.5 we know that as — oo 



sup (t + ^At) — /icq|L 

tG[0,T] 







and for any / G C{E) and i,j & Q 

sup \hf{t + tN){f,fif{t + t^)) - hf{t + tr,){f,fif{t + t^))| 

tG[0,T] 

Note that this also implies that if L is a function in the class T then 

sup |L(/i^(t + t7v))| ^ as iV ^ oo. 

te[o,T] 



0. 



(4.146) 



(4.147) 



(4.148) 



(4.149) 



We argued before that the sequence of processes {u^ : G N} is tight. Let {i^{t) : t > 0} 
be a limit point. Then along some sequence k^, ^ u as N oo. Since =^ oo, using 
the continuous mapping theorem and (4.147) we obtain that along the subsequence kN 



sup 

te[o,T] 



as A^ -> OO, 



(4.150) 
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where ( is any of the continuous functions (p, (p, Cf' ' Cf' 5 Cf' ' Cf' defined in Remark 4.9. 
From (4.142), (4.146), (4.149), (4.150) and (4.141) we get that along the subsequence kN 



sup 

te[o,T] 



/ Af F(/i^(s))rfs - / AiF{v{s))ds 
Jo Jo 



as iV ^ 00. (4.151) 



Using (4.151) and the continuous mapping theorem we can conclude that for any F G Co, as 
— )■ 00, the sequence of martingales rrip (given by (4.145)) converges in distribution along 
the subsequence /cat to the martingale given by 

F{u{t))-F{u{0))- [ AiF{u{s))ds. (4.152) 
Jo 

This shows that {z/(t) : t > 0} is a solution to the martingale problem for A^. Let vr G 
V{V{E)) be the distribution of r(//(0)). Since /i^(0) =^ /i(0) as — )■ 00 and F is a continuous 
map we must also have that z/^(0) =^ '^(O), where z/(0) has distribution it. We argued in 
Section 2.4 that the martingale problem for each A; is well-posed. Hence {z^(t) : t > 0} is 
the unique solution to the martingale problem for (A;,7r) and thus ^ u a.s N 00, 
along the entire sequence. Moreover the limiting process has sample paths in C-p(F)[0,oo) 
almost surely. 

Let {/i^(t) : t > 0} be the process defined by (2.43). Pick a i G Q and / G C{E). From 
(4.137) for any < t < o"^ — t^v "we can write 

</,/^fW> = (/,/if(i + ^^)> 

+ J2 (hfit + tN){f, /if (t + t^)) -h^{t + t^){f, /if (t + t^))) A,(/.^(t + t^)). 

Since =^ u, =^ 00 and 0, (4.147) , (4.148) and the continuity of the sample 

paths of {i^it) : t > 0} imply that for any T > 

sup |(/,/2f (t)) - /ieq,i(/,i/(t))| ^ as AT ^ 00. 

tG[0,T] 

This holds for any / G C{E) and i E Q. Hence by Theorem 3.7.1 in Dawson[7], /i^ =^ hcqiy 
as A^ —J- 00 in the Skorohod topology on Dj^i (^p^ [0, 00). This completes the proof of Theorem 
2.6. ^ □ 



Remark 4.10 In the statement of Theorem 2.6 we did not specify how the initial distribution 
7T of the limiting Fleming-Viot process {I'it) : t > 0} is related to the distribution of /i(0). 
However the above proof makes it clear that n is the distribution o/F(/i(0)) where F is the 
map defined by (4.119). 
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A Appendix. 

Lemma A.l For each h G let A{h) G M(g, g) be the matrix given by (2.13). Suppose 



that Assumption 2.1 is satisfied and let h, 



eq 



be the nonzero vector such that 



^(/ieq)/ieq = Oq 



(A.l) 



Then we have the following. 

(A) The vector h^q is in R^, that is, all its components are strictly positive. 

(B) The matrix A{hcq) has as an eigenvalue with algebraic multiplicity 1. All the other 
eigenvalues of A{hcq) have strictly negative real parts. 

(C) There exists a unique vector t>cq G such that Weq^(^eq) = Og and {vcq, /ieq) = 1- 



(D) Let G{heq) G M(g, g) be the matrix given by 



G(h 



eq) 



h 



cq 



^eq) 



A{K 



eq7 



Then the matrix G{hcq) has the same eigenvalues as the matrix y4(/ieq). 

(E) Let G{hcq) G M(g — 1, g — 1) be the matrix defined by 

Gij{Kq) = Gij{hcq) - Giq{hcq) for all i,je{l,...,q- 1}. 

Then the matrix G{heq) is stable, that is, all its eigenvalues have strictly negative real 
parts. 

Proof. Observe that all the non-diagonal entries of the matrix A{hcq) are nonnegative. Such 
matrices are sometimes referred to as Metzler-Leontief matrices in mathematical economics 
(see Section 2.3 in Seneta [32]). Their important property is that they can be transformed 
to a nonnegative matrix by adding a constant multiple of the identity matrix. This allows 
extensions of the Perron-Frobenius type results to such matrices. 

Note that the matrix A{heq) is irreducible (part (C) of Assumption 2.1) and has as 
an eigenvalue with h^q as a right eigenvector (see (A.l)). Theorem 2.6 in [32] proves parts 
(A),(B) and (C) of this lemma. 

We now prove parts (D) and (E). Let P and its inverse be the matrices in M(g, g) 
given by 



— T 

0,-1 1 



and P 



f 



-1 



9-1 



■T . 



Observe that 1 G(/i, 



cq; 



Oq and 



P^Gih,q){P 



G{hcq) V 
—T 

0,-1 







(A.2) 
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where v is some vector in R^^^. Let L = Diag(/ieq) and Q = p-^L-^. Note that Qh 



P ^hcq = P ^Iq = Oq and hence 



eq 



This shows that QG{hcq)Q^^ = QA{hcq)Q~^ . Hence the matrices G(/ieq) and A{hcq) are 
similar and have the same eigenvalues. Part (B) of this lemma and (A. 2) imply that the 
matrix G{hcq) is stable. This completes the proof of this lemma. □ 

Let A, 6, ip and t/gq be as defined in Section 4.3 (just prior to Proposition 4.7). 

Lemma A. 2 Fix a to > 0. For {h,t) G t/cq x [0,to] consider the following matrix equation 

<^{h,t,to) = Iq+ A^{4j{h,u))<l>{h,u,to)du, (A.3) 



where Ig is the q x q identity matrix. This equation has a unique solution in C'^[Ueq x 
[0, to]; MiK(g, q)) that satisfies the following for any h G f/eq, s > and < t < to- 

(A) Ifse[t,to] then^{h,t,to) = ^{h,t,s)^{h,s,to). 

(B) For any s > 0, $(^(/i, s),t,to) = $(/i,t + s,to + s). 

(C) Ifvo G Ml then <l>{h,t,to)vo G 

Proof. The function ip is in C^(M'^ x M_,_,M'^) and the matrix-valued function {h,t) t— t- 
A'^l^lh, t)) is in C^(M'' x]R+, M(g, q)). Standard existence and uniqueness results for ordinary 
differential equations guarantee that there is a unique solution for (A.3) in the class C^(f/eq x 
[0,to],MK(g,g)). 

Part (A) of the lemma is just the Chapman-Kolmogorov property (see Proposition 2.12 
in Chicone [G]). Note that due to the semigroup property for ip (similar to (2.19)) both 
^{2p{h, s),t,to) and ^{h,t + s,to + s) satisfy the same equation for t G [0,to]. Hence by 
uniqueness of solutions, part (B) is immediate. 

We now prove part (C). Note that only the diagonal elements of the matrix A^ {^{h,t)) 
can be negative. For t G [0,to] let 

c(t) = - min Aii(^(/i,t)). 

l<i<q 

Fix a Wo G and define 

"to 

c{s)ds ) ^{h, to — t, to)vo. 



L{t) = exp I / 

\Jto-t 



Then 

^ = (c(to - t)Iq + A^'i^ih, to - t))) L{t). 

But the matrix ^c(to — t)Iq + ^{'4^{h,to — t))j has all entries positive for any t G [0,to]- 

Since L(0) = vo G M^, for any t G [0,to] we have dL{t)/dt G M\. Therefore L{t) G and 
this proves part (C). □ 
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